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Abstract. We continue the study of the correlation functions for the point stochas- 
tic processes introduced in Part I (G. Olshanski). We find an integral representation 
of all the correlation functions and their explicit expression in terms of multivariate 
hypergeometric functions. Then we define a modification ("lifting") of the processes 
which results in a substantial simplification of the structure of the correlation func- 
tions. It turns out that the "lifted" correlation functions are given by a determinantal 
formula involving a kernel. The latter has the form (A(x)B(y) — B(x)A(y))/(x — y), 
where A and B are certain Whittaker functions. Such a form for correlation functions 
is well known in the random matrix theory and mathematical physics. Finally, we 
get some asymptotic formulas for the correlation functions which are employed in 
Part III (A. Borodin and G. Olshanski). 
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Introduction 

In this paper we compute the higher correlation functions of a remarkable family 
of stochastic point processes introduced by G.I Olshanski in [O]. The main result is 
an explicit integral representation for the correlation functions. We also introduce a 
modification of the processes that provides a connection with certain class of known 
processes. The present paper is the continuation of [O]. In the next paper [BO] we 
give some applications of our results. 

In general, our point processes have a representation-theoretic origin, see [KOV], 
[O]. Here, to be concrete, we give a shorter definition of the processes via finite-point 
approximations, see [O]. 

Let us consider the Thoma simplex Q (see [O]) 

0= | ai > a 2 > ... > 0, ft >(3 2 > ... >0; J] (a* + ft) < l} 
and its face 

fio = {«i>«2>..->0, ft>ft>...>0; J](ai + A) = l}- 

We define a sequence P ™)\ n = 1,2, . . ., of discrete probability measures on fi 
depending on two complex parameters z, z' as follows. The support of the nth 
measure consists of all points of the form 

f pi + 1/2 p d + 1/2 gi + 1/2 Q d + l/2 

\ n n n n 

d 

pi, gi,... p d , qd e {0,1,2,...}; d=l,2, ...; J^(pi + ft + 1) = n. 

i=i 

Here all ctj, for z > d are zeros. 
The measure of such a point equals 

Fzz ' {m ~ (t) n l\ fe!) 2 fe0 2 

where t = zz' and (a)k = T(a + k)/F(a) stands for the Pohgammer symbol. For 
the parameters z and z' there are two possibilities, if we assume that the measure 
of a point described above is always strictly positive, 

(1) z' = z and z <£ Z 

(2) z and z' are real and m < z, z' < m + 1 for some m G Z. 

One can prove that we really obtain a sequence of probability measures that 
weakly converges as n — > oo to some probability measure P 2Z / on O. 
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The measure P K Z J can be regarded as a probability distribution on the set of 
all Young diagrams of order n, if we consider the numbers (pi, . . . ,Pd] Qi, ■ ■ ■ ,Qd) 
as the Frobenius coordinates of a Young diagram, see [M, §1.1] for the Frobenius 
notation. Furthermore, the values P^) (p\q) can be regarded as 'moments' of the 
limit measure P zz > . Namely, 



where (p\q) is considered as the Frobenius notation for a Young diagram A; dim A 
is the dimension of the complex irreducible representation of the symmetric group 
S n corresponding to A (see [M], [JK]); and s\ stands for the so-called extended 
Schur function, see [O] , [KOO] for details. These moments lie in the base of all our 
computations. 

Now we shall introduce the point process, see [DVJ] for general information 
about point processes. Let us denote by I the punctured interval [0, 1] \ {0}. We 
shall call a set that does not intersect a sufficiently small interval (— e, e) a test 
set. A configuration is, by definition, a finite or countable system of points in I 
such that its intersection with any test set is finite. For any point of O we define 
the corresponding configuration (cki, ai2, • • • , —@2, • • • )> w h ere we suppose that 
ctj and Pi are nonzero. Then the probability measure P zz > on O provides us a 
probability measure on the configurations, i.e. a point process. We compute the 
correlation functions of these processes. 

In what follows we use the notion of the correlation measure as well as that of 
the correlation function. The latter is defined as the density of the former with 
respect to the Lebesgue measure. As a distribution it always exists. Note that in 
[DVJ] the correlation measure is referred to as 'factorial moment measure '. 

We shall work with distributions all over this paper. Everywhere below the word 
'distribution' is understood as a synonym of 'generalized function', and our notation 
for the distributions coincides with the usual one for generalized functions. 

The paper is organized as follows. 

Chapter 1 contains the computation of the correlation functions. In [O] the com- 
putation was reduced to a rather complicated multidimensional moment problem, 
which we formulate in Section 1.1. This moment problem contains the values of 
irreducible characters of finite symmetric groups. In Section 1.2 we give a refor- 
mulation of the well-known Murnaghan-Nakayama rule for computing these values. 
The results of Section 1.2 allow to reduce the moment problem to a finite num- 
ber of simpler moment problems; this reduction is carried out in Section 1.3. In 
Section 1.4 we introduce three simple operations on compactly supported distri- 
butions, these operations correspond to some natural operations on the moments 
of distributions. Using these operations, in Section 1.5 we construct the solutions 
of simplified moment problems stated in 1.3. In this construction we use certain 
distributions Id, d = 0, 1, 2, . . . with given moments. To clarify the idea of construc- 
tion in 1.5, we postpone the explicit computation of i^'s to Section 1.6. After this 
computation we already have an expression of the correlation functions in the form 
of a finite sum of some known distributions. However, the summation is taken over 
a complicated combinatorial set; and we substantially simplify the description of 
this set in Section 1.7. 

It turns out to be much easier to describe the correlation functions (and the 
whole process), if we restrict ourselves to the behaviour of a^s (or /%'s) only, see 




3 



the definition of the point processes above. We do this in Chapter 2. Section 2.1 
explains that in this case only one summand of the complicated finite sum from 
Chapter 1 is relevant. In Section 2.2 we produce a simple integral representation 
of this summand; and it awfully resembles that of multivariate hypergeometric 
functions. In Section 2.3 we introduce Lauricella hypergeometric functions of type 
B and in the next section we express the correlation functions via these Lauricella 
functions explicitly. These two sections are considered to be optional. They contain 
rather hard technical work, and results obtained there will be used in also optional 
Sections 4.2 and 4.3. In Section 2.5 we prove that our point processes are simple, 
i.e., that there are no repetitions of ct^'s or fa's with probability 1. 

Chapter 3 is devoted to a modification of our processes, which we call lifted pro- 
cesses. The idea is that a simple lifting of measures P zz > from Qq to the bigger 
space of pairs of convergent series with positive terms substantially simplifies the 
correlation functions. In Section 3.1 we introduce a general formalism of the lift- 
ing and show that its application to well-known Poisson-Dirichlet processes gives 
just a Poisson process. In Section 3.2 we demonstrate that the lifting reduces mul- 
tidimensional integral representations of correlation functions obtained before to 
determinants built from a function in 2 variables with a given integral represen- 
tation. In the context of Chapter 2 the situation is even more simpler: the lifted 
correlation functions are given by determinantal formulas with a kernel expressed 
via Whittaker functions, we prove this in Section 3.3. We call the kernel the Whit- 
taker kernel, the explicit formula for it is one of our main results. Determinantal 
formulas for correlation functions provides a connection of our processes with those 
arising in random matrices and with certain models of mathematical physics. We 
hope to clarify this connection in subsequent papers. 

In Chapter 4 we compute the asymptotics of both lifted and non-lifted correla- 
tion function at the origin, which happens to be the same. This computation can be 
considered as an application of our integral representations. But its main destina- 
tion is to provide a base for considering so-called tail processes, which describe the 
behaviour of particles 'infinitely close to zero', see [BO] for further explanations. It 
turns out to be much easier to compute the asymptotics in the lifted case: it suffices 
to work only with the Whittaker kernel. This is done in Section 4.1. Section 4.2 
and 4.3 are optional: they, again, contain many technicalities which do not produce 
surprising results. However, without them this work would be incomplete. In Sec- 
tion 4.2 we construct asymptotic formulas for Lauricella hypergeometric functions, 
and in Section 4.3 we use them to compute the asymptotics in the non-lifted case. 
The computation substantially exploits formulas from Section 2.4. 

This work was totally inspired by G. I. Olshanski. Besides that he set up the 
problem, numerous stimulating discussions with him played the crucial role for the 
whole work. I am happy to express him my deep gratitude. 

1. General Structure of The Correlation Functions 

1.1. The moment problem. As was shown in [O], the calculation of the cor- 
relation functions of our processes can be reduced to the calculation of the fam- 
ily of probability measures on cubes [— 1, l] n , n = 1, 2, ... , with their moments 
known. These measures will be referred to as controlling measures and denoted by 
a n 2 (dx\, . . . ,dx n ); we shall denote their densities with respect to the Lebesgue 

( zz ) 

measure by a n (xi, . . . ,x n ). The densities are considered as distributions (with 
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compact supports). Here z and z' are parameters of the process, see Introduc- 
tion. The moments of the controlling measures are given by the following explicit 
formulas, see [O, Proposition 3.3], 



/ 



r h ■ ■ ■r ln rT { - zz 'Udr, fir ) - \^ \^ (pi,~- ,Pd\qi,- ,Qd) 

x l x n°n l"^l> ••• i u ^n) ~ / J / , X(h+1,... ,l n + l) 



r— 1,11™ pi>...>p d >o 

<?!>.. .>q d >0 
\p\ + \q\+d=\l\+n 

(1.1) 

' ( Z+ l) pi (- Z+ l) pi ( Z ' + l) qi (- Z ' + l) qi ( 1 



X 



(*)|p| + M+dti \P* + ?7 + 1 

where h, . . . , l n = 0, 1, 2, . . . ; 

|p| =Pi + .. .+Pd, |<?| = 9i + • • • + qd, \l\ = h + • • • + l n + n; 

(pi, . . . , Pdl^i, • • • 7 %) is the Frobenius notation for a Young diagram, see [M, §1.1]; 
and Xmi,... ,m k i s an arbitrary Young diagram, mi, . . . , mk are positive integers) 
stands for the value of the irreducible character \ X of the symmetric group Sja| at 
any permutation with k cycles of length mi, . . . , m/-, see [JK], [M]. 

The connection between the correlation and controlling measures is provided by 
the following statement proved in [O, Proposition 4.3]. 

Theorem 1.1.1. On the set {{x\, . . . , x n ) £ I n \ Xi ^ Xj for i ^ j} 

Pn (-^l) • • • ) x n) = "j T°~n (^1) • • • ? x n)- 

\ x l x n | 

Moreover, using [O, Proposition 4.2] we shall demonstrate (see Section 2.5) that 
the set 

{(xi, . . . , x n ) G I n \ Xi = Xj for some i ^ j} 

is a null set with respect to the nth correlation measure. 

Thus, our problem is completely reduced to the calculation of the controlling 
measures. 

1.2. Combinatorics of the Murnaghan-Nakayama rule. In the following 
we shall need a reformulation of the well-known Murnaghan-Nakayama rule for 
computing the values of irreducible characters of finite symmetric groups. 

Theorem 1.2.1 (Murnaghan-Nakayama rule, see [M, §1.7, Ex.5]). Let A be 

a Young diagram, |A| = n; p = (pi, . . . , pk), all pi are positive integers, p\ + . . . pt = 
n. Then 

X X P = E(-!) ht(S) 

s 

where the sum is taken over all sequences of Young diagrams S = (X^ C A^ 1 ) C 
. . . C A( fc ) = A) such that A^ — A^ -1 ) is a skew pi-hook; ht(S) is the height of S, 

ht(S) = ^ht(A (i) - A (i_1) ). 

i 
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(The height of a skew hook is one less than the number of rows it occupies.) 
Let us introduce some notation. 

We shall consider certain positive integral-valued variables which will be called 
linear horizontal and linear vertical blocks. 

We also use variables whose possible values are pairs (p, q) G {0, 1, ... } 2 , and we 
call these variables hook blocks. 

The filling number of a linear block is its value. For a hook block we shall call p 
the filling p-number and q the filling q-number of this block. 

The value of a linear block will be also called the cardinality of the block. For a 
hook block we define its cardinality as the sum of its two filling numbers plus one. 

A fragment is, by definition, a partially ordered finite set of blocks subject to the 
conditions 

(1) there is exactly one hook block that precedes all the others; 

(2) all linear horizontal blocks are linearly ordered; 

(3) all linear vertical blocks are linearly ordered; 

(4) there are no other order relations. 

Finally, we shall call a finite set of fragments with linear ordering of all their 
blocks (the ordering is compatible with the partial orderings inside the fragments) 
a structure. 

A structure is called filled if all filling numbers of its blocks are known. 

The ordered set of cardinalities of all blocks of a filled structure is called the 
cardinality of the filled structure (the order of blocks' cardinalities is the same as 
the order of blocks in the structure) . 

We shall denote by v(T) the number of linear vertical blocks in T. 

Let us define the filling p-number (q-number) of a fragment as the sum of the 
p-number of its hook block and the filling numbers of its linear horizontal blocks 
(the sum of the q-number of the hook block and the filling numbers of its linear 
vertical blocks, respectively). 

Thus, any filled structure T with d fragments produces a set of pairs 

{(Pi,Qi),...,(P d ,Qd)} 

which consists of the filling p- and q-numbers of the fragments. 
Let us define the sign of T as follows 

sgnT = sgn(lJ(P, - P 3 ){Q t - Q s )) (-1) E Qj+v(T) . 

Our main statement in this section is the following 
Theorem 1.2.1'. Under the conditions of Theorem 1.2.1 

(1-2) X J = $>S nT ' 

T 

where the sum is taken over all filled structures with the cardinality p = (pi, . . . , pt) 
such that the sets {Pi, . . . , Pa} and {Qi, . . . , Qd} of the filling p- and q-numbers of 
the fragments coincide, up to a permutation, with the Frobenius coordinates of the 
Young diagram X. 
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Remark 1.2.2. For the calculation of the nth correlation function we shall consider 
the structures with n blocks. For n = 1 there is only one such structure; it has one 
hook block. For n = 2 we have already 3 such structures: the structure with two 
fragments, each of them consists of one (hook) block; the structure with one hook 
and one vertical block; and the structure with one hook and one horizontal block. 

Proof of Theorem 1.2.1'. We shall call the structures, which can be obtained from 
a given one by throwing out several last (in the sense of the ordering inside the 
structure) blocks, the substructures of the given structure. A filling of the initial 
structure induces the fillings of all its substructures. We shall call a filled structure 
proper if for all its substructures (including itself) the following assertion holds: all 
filling p-numbers (as well as q-numbers) of the fragments are mutually distinct. 
The statement of Theorem 1.2.1' immediately follows from the following lemmas. 

Lemma 1.2.3. All proper filled structures T which participate in the sum (1.2) are 
in one-to-one correspondence with the sequences S from Theorem 1.2.1. Moreover, 

(1.3) sgnT = (-l) ht(S) 

for corresponding T and S. 

Lemma 1.2.4. The part of the sum (1.2) corresponding to all improper filled struc- 
tures T vanishes. 

Clearly, these lemmas and Theorem 1.2.1 imply Theorem 1.2.1'. 

Proof of Lemma 1.2.3. Any removal of a skew 9- hook from a Young diagram is 
equivalent either to a removal of two Frobenius coordinates of the diagram with 
their sum equal to 9— 1, or to a reduction of a Frobenius p-coordinate or (/-coordinate 
hj9. 

More precisely, if we remove a #-hook s that does not intersect the diagonal and 
lies to the right of it occupying the rows from ith to jth, % < j, then it is equivalent 
to the following changes in the sequence of Frobenius p-coordinates {p\ > . . . > pa) 
of our diagram 

(p 1 > . . . > p d ) -> (pi > . . . > pi- t > p i+1 > . . . > pj > pi - 9 > p j+1 > . . . > p d ). 
By definition, 

(_l)ht(s) = 

If we remove a 6>-hook s that lies completely to the left of the diagonal occupying 
the columns from i'th to j'th, then, by analogy, it is equivalent to the following 
changes in the sequence of Frobenius g-coordinates of the diagram 

(q 1 > . . . > q d ) -> (q 1 > . . . > q v _x > q v+1 > . . . > q r > q v -9 > q r+1 > . . . > q d ) 

The cardinality 9 of s is the sum of the number of columns occupied by s plus the 
number of rows occupied by s minus one. Therefore, 

(_ 1 )ht( s ) _ 

Finally, if we remove a #-hook that intersects the diagonal of the diagram, oc- 
cupies rows starting from the ith one, and columns starting from the i'th one; it 
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is equivalent to removing the ith Frobenius p-coordinate and the z'th Frobenius 
g-coordinate, and 

Pi + qi' = - 1. 

This #-hook can be divided into three parts: p^-hook lying to the right of the 
diagonal, gi'-hook lying to the left of the diagonal, and a box lying on the diagonal. 
Using this and previous considerations, one can easily compute that 

(_l) ht (s) _ (_iy'-i+Qi' _ 

We define the correspondence from the hypothesis of Lemma 1.2.3 in such way 
that three operations described above correspond to deleting a linear horizontal, a 
linear vertical, or a hook block of a structure respectively. 

Formally it can be made by induction on k (the number of blocks in T or the 
length of p). Namely, for k = 1 the sequence = A^ ) C A* 1 ) = A, A = (p\q), 
corresponds to the structure T with one hook block, and its p- and q-filling numbers 
are exactly p and q. Furthermore, if a sequence S' with k = m is obtained from 
some sequence S with k = m+1 by throwing out the last skew hook \( rn + 1 ') — \( m \ 
then the corresponding to 5" structure T' is obtained from corresponding to £ 
structure T by throwing out the (m + l)th block. If the hook thrown out lies to the 
right of the diagonal, to the left of the diagonal, or intersects the diagonal, then the 
corresponding block is linear horizontal, linear vertical, or hook respectively. The 
Frobenius p- and (/-coordinates of A = A^ m+1 ^ coincide (up to permutations) with 
p- and q-filling numbers of fragments of T. 

Let us check the coincidence of signs (1.3). We shall do this by induction on the 
number of blocks. If our structure has only one (hook) block with p and q as its 
filling numbers then the corresponding diagram A = (p\q) is a skew hook itself, and 
the corresponding sequence S is = A^ ) C A^ 1 ) = A. In this case 

sgnT = (-l) ht(S) = (-1) 9 . 

Let us consider a structure T with m + 1 blocks such that the (m + l)th block of 
T is linear horizontal and its filling number is 9. Let us denote by T' the structure 
obtained by removing the (m + l)th block from T (as above). At the beginning of 
the proof we considered the corresponding removal of a skew #-hook s. The p-filling 
numbers (Pi, . . . , Pa) of fragments of T are permuted numbers {p\ > . . . > pa)- It 
is easy to see, that if we reduce one of P's corresponding to pi by 9, then the sign 
of n^jl-^i ~~ Pj) an d sgnT will be multiplied by (— l) 1- - 7 , that coincides with the 
sign (— l) ht ( s ). So we have reduced the coincidence of signs for T and corresponding 
sequence S to that for T' and corresponding sequence 5", and the number of blocks 
in T' is less then that in T. One can easily check that the same reduction works for 
the last block of T being linear vertical or hook as well. Modulo this check, (1.3) 
is proved by induction. 

The condition that all structures must be proper is necessary because the Frobe- 
nius p- or (/-coordinates of any Young diagram are pairwise distinct. □ 

Proof of Lemma 1.2.4- Let us exclude from our consideration improper filled struc- 
tures with equal p- or q-numbers of fragments; their signs are zeros. All the other 
improper structures will be divided into pairs with opposite signs. 
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Let us fix one such structure. Let us choose from its substructures the biggest 
with equal p- or q-numbers. To be concrete, let there be two equal filling p-numbers 
of some two fragments. Let us exchange in this fragments all linear horizontal 
blocks which distinguish the initial structure and the substructure. The set of such 
blocks is not empty, because the initial structure has different p-numbers. Thus we 
obtained a new filled structure which will be the pair for the initial one. Obviously, 
the procedure described above is involute, and the structures in one pair have 
opposite signs (the Vandermonde determinant Yl (Pi ~ Pj) ( or Yl (Qi ~Qj) m case 

i<j i<j 

of q-numbers) changes its sign). □ 

Remark 1.2.5. The coincidence of signs (1.3) for a structure without linear blocks 
was communicated to me by G. I. Olshanski. 

1.3. A simplification of the moment problem. Our goal in this section is 
to represent the nth correlation measure as a sum of distributions with compact 
supports; the summands will be naturally parametrized by (not filled) structures 
with n blocks. The distributions will be defined as the solutions of moment problems 
which are much simpler than the initial one, see 1.1. 

Proposition 1.3.1. The moments of the nth controlling measure have the form 

T 

U(z + l) Pi (z' + l) Pt (-z + l) Qi (-z' + l) Qt 



x — — — det 



1 



(t) E{Pi+Qi+1) nPiW yPi+Qj + i 

i 

where the summation is taken over all filled structures T of the cardinality {l\ + 
1, . . . , l n + 1); (Pi, . . . , Pd', Qi, ■ ■ ■ , Qd) o,re the filling numbers of the fragments of 
T , d is the number of the fragments. 

Proof. This statement is an immediate corollary of the formula (1.1) and Theorem 
1.2.1'. The only thing to mention here is that due to the well known formula for 
the Cauchy determinant (see [W]) 



det 



Yl{(Pi-Pj)(Qi-Qj)) 

i<3 



Pi+Qj + ij noPi+Qj + i) 



we may use the equality 



3 S ndet ( p. + n. + i ) = sgnn((^ - Pj)(Qi ~ Qi))- D 



Let us represent the sum from Proposition 1.3.1 as the double sum: the outer 
summation will be taken over all structures with n blocks, and the inner summation 
will be taken over all possible fillings of the structure given by the outer sum. Thus 
we obtain partitions of all the moments into summands parametrized by the unfilled 
structures with n blocks. In Section 1.5 we shall construct the distribution with 
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a compact support with its moments equal to the summands corresponding to a 

fixed structure T; this distribution will be denoted by a n (T) . Now we are going 

to reformulate the moment problem for o~n (T). 

Let us fix a structure T with n blocks and d fragments. Let us enumerate the 
fragments arbitrarily by the numbers from 1 to d. We denote the number of linear 
horizontal blocks in the ith fragment by \ii and the number of linear vertical blocks 
-by Vi . Thus 

d 

^2(pt + Vi + 1) = n 

i=i 

(in each fragment there is exactly one hook block). 
In this notation v(T) = Yli=i v i- 

We say that the variable xj corresponds to a block if the block has the number 
j in the total blocks' ordering of T. Let us rename the variables x\, . . . , x n by the 
letters 

fli, {HVU, W,Y'i ■ i = l,...,d; 

so that the variable corresponds to the hook block of the ith fragment; {bl}^ =1 

correspond to the horizontal blocks of the ith fragment; and {c J i }^L 1 correspond to 
the vertical blocks of the ith fragment. 

By analogy with this, let us rename the numbers (Zi,... , l n ) by the letters 
({^},{^},{C/}). 

Proposition 1.3.2. With the preceding notation 

T 

the sum is taken over all structures T with n blocks and 

i=i 

(1.4) =t d (-l) EQl+u(T) •- 



(t) E (p i+Qi +i)nPiW 

i i 
i=l,... ,d 

1 



x det 
where 



P l + Qj + l 



(1.5) 



Pi=A' i+ f:(Bi+i) 



3 = 1 



Proof. The only freedom in filling a structure T with a fixed cardinality is that we 
can vary filling p- and q-numbers of every hook block with their sum fixed. If we 
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denote the p-number of the hook block in the ith fragment of T by A\ and the 
q- number of the same block by A", then we observe that the cardinality of the 
block in question equals Ai + 1 = A\ + A" + 1. Computing the filling numbers of 

( zz') 

the fragments (denoted by Pi and Qi) together with the definition of ak (T) and 
Proposition 1.3.1 complete the proof. □ 

1.4. Contraction, pseudoconvolution, and diagonalization. In this section 
we shall introduce three operations on distributions which will be used further for 
the construction of solutions of our moment problems (see previous section for the 
setting of the problems) . 

Proposition 1.4.1 ('Contraction'). Let £2, ?7i, • • • ,Vk), k>0, be a distri- 
bution with compact support and the moments 

Then there exists a distribution g(£,rji,... ,rjk) with a compact support and the 
moments 

(g,Cvi (3l ---Vk f3k )= £ (-l) a2 m aia2/3l ...^. 

If, moreover, supp C {6 > 0, £2 > 0}, then in the domain {£ 7^ 0} the following 
equality holds 



f Ifl 

J \€\ + y 



Vk))dy. 



y>0 

We shall use the notation 



g(£, 771, ...,%) = C| i6 [(/ (£1, £2, m, • ■ • » Vk)] ■ 

Proof. The value of the distribution g on any test function r/i, . . . , ?7fc) is given 
by the formula 

£i^(£i,vi, -..,%) + 6^(-6, 7/1, . . . , 7/fc) . 



M) = {f, 



6 + 6 



Proposition 1.4.2 ('Pseudoconvolution'), see [O]. 

Let /(£i> • • • ? Cm) an d g(vi, • • • 7 ?7m) distributions with compact supports and 
the moments 

(/, 6" 1 • • • , £m° m > = m ai ... 0m , (<7, r// 1 • • • , ijj 3 ") = nft...^. 

T/ien i/iere exists a distribution h((i, . . . , Cm) compact support and the mo- 
ments 

(h, Cl 71 5 Cm 7 '™) = m 7i...7 m n 7i...7 m - 

Moreover, 

supp ft C supp / • supp 0. 

VFe s/iaW use the notation 

h = f&g. 

Proof. The value of the distribution h on a test function ^ is given by the formula 

„ m 

= J /(6,--- ,Cm)^(^7i,--- ,Vm)ip(£,iVi,--- ,£,mVm)Yl d & dr li- 



i=l 
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Proposition 1.4.3 ('Diagonalization'). Let /(£,?7i, . . . , rjk), k > 0, be a distri- 
bution with the moments 

Then the moments of the distribution 

g(£i, 6, m, • • • , Vk) = m, • • • , Vk)K& - Ci) 

We shall use the notation 
Proof. Obvious. 

(zz') 

1.5. The distribut ions o~ n (T) . In the present section we shall construct the 
solutions of the moment problems (1.4) for all T starting from a certain family 
of distributions. We shall also establish that only few of these solutions really 
contribute to the correlation functions. 

The crucial role in the whole construction is played by the family of distributions 
called &'{'■> ■ ■ ■ ', a d ), d = 1,2, ... , whose moments, by definition, are equal 

to 



1 



(I d (a[, a'{; ...;a' d , a» d ), (a'^ (a'/) A " • • • {a' d ) A \a'^) 
U(z + 1)a>(z' + Va'S-z + l) K {-z' + 1)a» 

(L6) = ~ rppfK^^j^ de \ A i+ A "+ 1 )' 

i 

The following statement will be proved in the next section, and now, in order to 
clarify the construction, we shall assume that it is true. 

Proposition 1.5.1. There exist distributions Id(a>i, a'{; . . . ;aa, a d ), d = 1,2,..., 
such that their moments are given by (1.6). Moreover, 

d 

(1.7) supp I d C {K, a'{; ...;a' d ; a'fi G R 2 + d \ + <) < !}• 

i=i 

We use the notation introduced in 1.3. Let us denote (for alH = 1, . . . , d) 

n po,i 



(i.8) ^=rKp' ^=rKr c ^ 



Moreover, let us denote by M.{ (for all i = 1, . . . , d) the operation of multiplication 
of a distribution by (e^) M< {a>'l) Vi ■ 
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Proposition 1.5.2. For any fixed structure T the compactly supported distribution 



o~n Z (T) m variables 



defined by 



^ z '\T) = t d (\{c^ (n^j ^n^*) 

gives a solution of the 'simplified moment problem 1 (1-4) ■ 

Proof. The following proof is nothing but sequential applications of statements from 
the previous section. Let us introduce the distributions 



X 



Y[M t ■I d (a' lt a'{;... ;a' dt a'£) 

\i=l J 



in variables aj, a", i = 1, . . . , d; 

in variables a[, a'-, {6- j^, {4}jLi, i = 1, • • • , d; and 

z=(n Ci ).K 

in variables c^, {fr^}^, {cj}y =1 , i = 1, . . . , d. Let us introduce the following nota- 
tion. For a distribution ... , we shall denote by m/ 1; . mmj i k (/) its (/i, . . . , 
moment. That is, 

mi ll ..., lfc (/) = (/,eJ 1 ...^>. 

Denote 

^' = {4}?=i> ^' = {^}ti, A = {A i }f =1 ; 

b = {b{}1% =1 , C = {Cj}tl j=1 ; 
v = {M?=i> f = {^}ti; 

cf. 1.3. Then, by Proposition 1.4.1, 

(1-9) m A , B ,c( z )= E (-1) |A " Ito ^'^",b,c(^) 

A'+A"=A 
A'i,A">0 

where the absolute value sign for a vector stands for the sum of its coordinates. 
Furthermore, by Proposition 4.3 we get 

(1-10) mA>,A»,B, C {Y) = (-^ l ^A' + B,A"+c( X ) 
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where 

Note that the sign (— l)\ c \ appears by the following fact: for a one dimensional 
distribution g(rf) 

(g(v),v l ) = (-i) l (g(-v),v 1 )- 

In our case we apply this fact to all variables {c^}, see the minus sign before cf in 
(1.8). 
Finally, 

m A>+B,A"+c(X) = m A , + B +fJijA „ + c +u (Id) ■ 

Recall the notation (1.5). It can be reformulated as 

P = A' + B + n, 
Q = A" + C + ls. 

In particular, (-1)I A "I + I^| = (_l)IQI+M = (_i)IQI+«(t)_ Combining (1.9), (1.10), 
(1.11) together we obtain (using (1.6)) 

mA,B,c(Z)= £ (-l) IQI+w(T W**) = £ (_1)IQ ' 



. . (t)\P\ + \Q\+d 
A'+A"=A A'+A"=A 



A'i,A'/>0 



* (_z+ l) Pi {z> + l) Pi (-z + l) Qt {-z' + \) Qi 
l\ P^l ' de \p i + Q j + l 

Hence, the moments of t d Z coincide with those from (1.4), as was to be proved. □ 

Let us call a structure substantial if there is at most one linear block in each 
its fragment. In the preceding notation it means that \ii + vi < 1 for all i. The 
next statement shows (see also Theorem 1.1.1 and a few words after it) that only 
substantial structures are relevant. 

Proposition 1.5.3. If a structure T is not substantial then 

suppa { n zz '\T) c {(zi,... ^in^-n^ - ^) =°}- 

i i<j 

Lemma 1.5.4. Let the (a±, . . . ,a.k)-moment of a distribution . . . , £&) with 
compact support depend on a±, ct2 only via their sum. Then 



supp/ C {(&,.. . ,Cn)| 6=6}- 
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Proof of Lemma 1.5.4. (fi - i, • • • , ffc) = 0. 

Proof of Proposition 1.5.3. Let us use the preceding notation. If a structure T is 
not substantial then there are three possible cases: Hi > 2 for some i, Vi > 2 for 
some i, or there exists such % that ^ = i/j = 1. 

Let us first assume that d = 1, i.e, our structure consists of only one fragment. 
We shall omit the subscript % of all variables and operations introduced above, 
because in this case % is identically equal to one. 

Suppose /j, > 2. Then the moments of \t), see (1.4), depend on B 1 and B 2 
only via their sum (because they depend on P = B 1 + B 2 + . . . ). Thus, Lemma 
1.5.4 implies the assertion of the theorem. 

Suppose v > 2. Then, again, the moments of a^ z \t), see (1.4), depend on C 1 
and C 2 only via their sum (because they depend on Q = C 1 + C 2 + . . . ) . Again 
Lemma 1.5.4 completes the proof. 

Finally, suppose /j, = v = 1. By Proposition 1.5.2 we know that 

Mh(a, a") =X, VX = V D" X = Y, CY = Z, 

a^'\T) = tZ. 

We have 

X = Mh(a', a") = a'a"h(a', a"), V" X = a'a"h(a', a")S(c + a"), 

Y = V'D"X = a'a"h(a', a")6(c + a")8(b - a'). 

Furthermore, applying Proposition 1.4.1 (we can use the explicit formula for the 
contraction because the support of Ii(a',a") lies in the domain a', a" > 0, see 
Proposition 1.5.1), for a ^ we obtain 



Z(a, b,c) = C- Y(a', a", b,c) = [ (Y(a, y, b, c) + Y(y, -a, b, c))dy 

Jy>0 \ a \ + V 

= I | r~ (ayl 1 (a,y)6(c + y)6(b-a)-ayl 1 (y,-a)6(c-a)6(b-y))dy 

Jy>0 \ a \ + V 

tt 2 c ci 2 b 

ii(a, — c)S(b — a) + Ii(b, —a)5(c — a). 



a — c b — a 

Obviously, the support of the last expression lies in the domain 

{(a, Ml (6-a)(c-a) = 0}, 

as was to be proved. 

In case of arbitrary d we use the fact that the operations Ci, A4i with different 
subscripts act on different sets of variables. That is why the same considerations as 
above applied to a fragment with linear vertical and linear horizontal blocks prove 
the assertion. □ 

Now let us fix a numeration of the fragments of a substantial structure T (this 
numeration was arbitrary up to this moment) in such a way that 



A*< = 1, v% = 
Hi = 0, Vi = \ 
Hi = 0, v t = 



i = l,...,m 1 , 

i = mi + 1, . . . , mi + m 2 , 

i = mi + m2 + 1, • • • , mi + rri2 + = d. 
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In other words, we denoted the number of the fragments with horizontal blocks by 
mi and put them first, then we put fragments with vertical blocks, and, finally, 
we have 7773 fragments without linear blocks at all. 
According to Theorem 1.1.1 we also set 

(1-12) P { : z '\T) = -A ZZ 'XT). 

\ x l x n I 

It turns out that the case 7773 = is much simpler than the general situation. 
Namely, in this case the formulas for On* ^ (T) and for pn Z (T) are rather simple. 

Theorem 1.5.5. Let 7773 = 0. Then in the domain where all the variables are 
nonzero and pairwise distinct 

n 

(1-13) (T) = t d J] — r-r/dCri, -s i; . . . , r d , s d ) 

I i ~T \i>i\ 

where 

(ri,... ,r d ) = (61,... , 6 mi , a mi+ i, . . . ,a d ); 
(si) • • • ,Sd) = (ai, • • • , a mi , c mi+ i, . . . , Cd). 

Remark 1.5.6. Theorem 1.5.5 and Proposition 1.5.1 imply that for 7773 = 

Pn (T) is a distribution defined in the domain where the sums of the absolute 
values of variables corresponding to the two blocks of each fragment are nonzero 
(each fragment contains exactly two blocks). Moreover, Theorem 1.5.5 implies that 

(zz r ) 

the support of p y n ' (T) lies (note that we have thrown out the set Yli<j ( x i~ x j) = 0) 
in the set where all variables corresponding to linear horizontal blocks (i.e. {bi}) 
are nonnegative, all variables corresponding to linear vertical blocks (i.e. {q}) are 
nonpositive, and every two variables corresponding to the two blocks of the same 
fragment (i.e. and bi for % < 777-1, or eij and q for m\ < i < mi + mi) have 
different signs. 

Proof of Theorem 1.5.5. We shall work out a detailed proof in two cases: (mi = 
1, 7772 = 0), and (mi = 0, 7772 = 1), and then give an explanation how to get the 
proof in the general case. 

Suppose 777-1 = 1 and 7772 = 0. We shall follow all the steps demonstrated in the 
proof of Proposition 1.5.2. By this proof we know that (cf. the proof of Proposition 
1.5.3) 

Mh(a', a") =X, VX = V'X = Y, CY = Z, 
ai zz '\T)=tZ. 

We have 

X = Mh(a', a") = a'hia', a"), Y = V'X = a' h(a' , a")S(b - a'). 
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Using Proposition 1.4.1, namely, the explicit formula for the contraction, for a ^ 
we get 



Z(a, b) = C a a , a „ ■ Y(a', a", b) = [ -J^-(Y(a, y, b) + Y(y, -a, b))dy 

Jy>0 \ a \ + V 

■ (aii (a, y)6(b - a) + yh(y, -a)6(b - y))dy 

= [ — t — h(a,y)dy ■ 8(b- a) - -^—h(b,-a). 
J y >o a + V b-a 





a 




\a\ 


+ y 



Note that the support of the first summand lies in the set where a = b. Thus, we 
have showed that in the domain where a ^ and a ^ b 

p[ zz ' \T) = tZ(a, b) = t-^-h(b, -a), 



as was to be proved (in this case r\ = b and s\ = a). 
Suppose mi = and 1112 = 1. Then 



Mh (a, a") =X, VX = V" X = Y, CY = Z, a ( 2 zz ' } (T) = tZ; 
and 

X = Mh(a', a") = a"h(a', a"), Y = V" X = a"h(a', a")5(c + a"). 
Again, by Proposition 1.4.1, when a^Owe get 

Z(a, c) = O a „ ■ Y(a', a", c) = [ -^—(Y(a, y, c) + Y(y, -a, c))dy 

jy>o \ a \ + y 

= I i FT {yh(a,y)8(c+y) - ah(y,-a)6(c- a))dy 

Jy>o \ a \ + y 

CLC f oF 1 

= h(a,-c)+ h(y,-a)dy ■ 5(c- a). 

a - c J y > y-a 

Now the support of the second summand lies in the set where c = a, and hence in 
the domain where a ^ and a^cwe have 

p { 2 zz '\T) = tZ(a, c) = t-^-h(a, -c). 



In this case r\ = a and s\ = c. 

In the general case for arbitrary number of fragments d we can apply the con- 
siderations demonstrated above to each of the fragments. Indeed, the operators 
Ci, T>i, M.i act on different variables if they have different subscripts. Moreover, 
each of the fragments is of one of the types considered above: it either contains a 
linear horizontal block (like in the case m\ = 1, m>2 = 0) or a linear vertical block 
(like in the case mi =0, m 2 = 1). Thus, in general case the proof is obtained by 
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word for word applications of one of our two previous considerations to the appro- 
priate fragments: the first part applies to the first mi fragments, and the second 
one - to the last m,2 fragments. □ 

The case m.3 > is a bit more complicated. As we have seen in the proofs 
of Proposition 1.5.3 and Theorem 1.5.5, the considerations for different fragments 
are independent. That is why, in order to understand the general situation, let us 
consider the unique structure that consists of only one (hook) block. That is, in 
our notation, mi = = 0, = 1. Then 

Mh(a', a") = I (a', a") = X, VX = X = h(a', a") = Y, 

CY = Z, a[ zz '\T) = tZ. 
By Proposition 1.4.1, for a ^ 

Z{a) = C a a , a „ ■ Y(a\ a") 



and in the same domain 
(1.14 p\ '(a)=t ■ dy + t — — dy. 

J y >o a + y J y >o m + y 

(Note that we have completely calculated the first correlation function, because 
there exists only one structure with one block.) Proposition 1.5.1 implies that the 
supports of the two summands do not overlap: in the first summand a > 0, while 
in the second a < 0. 

The answer (1.14) can be reformulated in the following way. Consider all sub- 
stantial structures with 1 fragment, which contain the initial structure as a sub- 
structure, and such that for them ms = 0, i.e., they have a linear (vertical or 
horizontal) block. In our case there are exactly two such structures - those that 
were considered in the proof of Theorem 1.5.5. Let us fix one of these two structures 
and denote it by T. Then it contains one extra block with respect to the initial 

(zz') ^ 

structure T. We take the distribution p 2 (T) described by Theorem 1.5.5 and 
integrate it over the variable corresponding to the extra linear block. Let us denote 
the answer by p\ zz (T, T). Then we have proved (formula (1.14)) that 

(1.15) p[ ZZ '\T) = J^p[ zz '\T,f) 

f 

where the sum is taken over all 'enveloping' structure described above (i.e., the 
sum contains two summands). As was mentioned before, the supports of the two 
summands do not overlap. 

This complicated explanation of (1.14) has only one advantage: an analog of 
formula (1.15) holds in general case. 

Let T be an arbitrary substantial structure. We fix the numeration of its frag- 
ments as described before Theorem 1.5.5. 
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a 



y>0 H +y 



(Y(a, y) + Y(y,-a))dy 



= I i \ a } (h{a,y) + h(y,-a))dy, 
J y >o m+y 



We shall call T an enveloping structure of T if T is substantial; T is a substructure 
of T; T and T have the same number of fragments; the numeration of the extra 
blocks in T (with respect to T) is compatible with the numeration of fragments of 
T; in each fragment of T there is exactly one linear block. 

In other words, in order to obtain an enveloping structure of T we have to add 
in each fragment of T with only one (hook) block a linear (horizontal or vertical) 
block, observing the ordering of fragments. The number of enveloping structures 
equals, obviously, 2 ma where, as before, 1113 is the number of fragments of T without 
linear blocks. 

For T being an enveloping structure of T we define 

(1.16) p^'\T,f)= J pi z *l s (f)dx n+1 ---dx n+m3 . 

2-n + li--- ; 2- ?i + m 3 

I.e., pn Z \t,T) is obtained from Pn+m 3 (T) by integrating the latter over all 

variables which correspond to the extra ms linear blocks added to T for obtaining T, 

/ t\ 

and Pn+rn 3 (T) is given by the general formula (1.13) (by definition, any enveloping 
structure has no fragments with only one block). 

(zz') ^ 

Remark 1.5.7. The distribution p n (T,T) does not depend on the numeration 
of fragments in T fixed before. A change of this numeration is equivalent to a 
change of the numeration of extra variables (x n +i, . . . ,x n + m3 ) that does not 
affect the result of the integration in (1.16). 

Theorem 1.5.8. In the domain where all variables are nonzero and pairwise dis- 
tinct 

(1-17) ti z '\T) = Y,P { : z '\T,T) 

f 

where the sum is taken over all enveloping structures T of T . 

Remark 1.5.9. Using Remark 1.5.6 and (1.16) it is easy to see that the supports 
of different summands of (1.17) do not overlap (in the domain where all variables 
are nonzero and mutually distinct), the support of each summand lies in one of 
the hyperoctants {xi > or < 0, % = 1, . . . , n}, and the hyperoctants are different 

for different summands. Thus, the support of pn (T) lies in the union of 2 7713 
hyperoctants which are described by the following conditions (cf. Remark 1.5.6): 
all variables corresponding to linear horizontal blocks are positive, all variables 
corresponding to linear vertical blocks are negative, and every two variables corre- 
sponding to two blocks of the same fragment have different signs. 

Proof of Theorem 1.5.8. We have already done all the work. One have to apply the 
considerations, which led to the formulas (1.14) and (1.15), to each of fragments 
without linear blocks. The other mi + fragments are considered as in Theorem 
1.5.5. As we know, different fragments are treated independently, and the result of 
all these applications will be exactly the formula (1.17). □ 

In Theorems 1.5.8 and 1.5.5 we have completely determined (modulo the con- 
struction of Id, see next section) the correlation functions in the domain where all 
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variables are nonzero and mutually distinct, using the notion of structures. But the 
combinatorics of substantial structures and their enveloping structures seems rather 
complicated. In Section 1.7 we shall interpret the pairs (T, T) as the elements of a 
simpler combinatorial object. 

1.6. Construction of I d . The goal of the present section is to prove Proposition 
1.5.1 and to provide explicit formulas for the distributions Id, d = 1,2, .... We 
shall start using this formulas in Chapter 2. 
Let us introduce the distributions 

u a 

<t> a {u) = - + , ueR, aeC. 

1 (a + 1) 

For 3?a > —1, <p a , by definition, equals u a /T(a + 1) for u > and vanishes for 
u < 0, so it is a locally integrable function. For < —1, 4> a is defined via analytic 
continuation. For example, 0_i(tt) = S(u). 
We shall also deal with products of the type 

(f>ab = (f> a {u)(f>b{l - it), a,beC, 

which are well defined for the reason that possible singularities of the factors do 
not overlap. 

Let us also remind that in Proposition 1.4.2 we have defined an operation 
on compactly supported distributions in the same number of variables and called 
it pseudoconvolution. The characteristic property of the pseudoconvolution of two 
distributions is that its moments equal the products of moments of initial distribu- 
tions. 

Proposition 1.6.1. For every d = 1,2,... there exists a compactly supported 
distribution Jd(a>i, a>'{', • • • ; a' d , a d ) with moments 

(J d (a[, a'!; . . . ; a' d , a'j), (a[) A ^ (a'l) A " • • • W d ) A ' d W) A *) 

U(z + 1)a[(z' + 1)a[(-z + 1)a[(-z' + 1) A '/ d x 

(L18) = ~ n^!^!WE(^ + A» + i) U A i + ^ +1 ' 

i 

Moreover, 



(1.19) supp J d C {(a' l7 a'{; ...;a' d ; a'fi e R% d \ J^M + <) < 1 }- 

i=i 

Proof. Using pseudoconvolution, let us define the distributions Jd as follows 

(1.20) J d = F 1 QF 2 Q F 3 , 
where 

d d 

FiK, a'l- ...;a' d , a'fi = T(t) J] ^{a'M-AO • <Pt-d-i (l - J2« + a ")Y> 

i=i i=i 
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F 2 (a' 1 , a"; . . . ;a dl a d ) =Y[<i>z,-z-i(a'J<l>- ZtZ - 1 (a"); 

i=l 

d 

F 3 (a[, a'(; . . . ; a' d , a d ) = S(a[ - a'{, . . . , a' d - a d ) JJ X[o,i] ( a D 



i=i 

Here 

X[o,i]{u) = 



1, «6 [0, 1] 
0, u i [0, 1] 



is the characteristic function of the segment [0,1]. By the fact that the support 
of pseudoconvolution is a subset of the pointwise product of the supports of the 
factors, the inclusion (1.19) follows from the following obvious relations 

d 

supp F l C {K, a'(; . . . ; a' d ; a'£) G R™\ ^« + <) < 1}. 

i=i 

suppF 2 , suppF 3 c {(ai, a"; . . . ; «4 a' d )\ < a-, a-' < 1, i = 1, . . . , d}. 

Let us compute the moments of J a and compare them with (1.18). We shall 
use the following lemma which is a well-known generalization of the Euler's beta- 
integral (the integral below is also called Dirichlet integral). 

Lemma 1.6.2. For any a Q , a\, . . . , a m e C 

i 

0a o ( u o) ■ ■ ■ 4>a m (u m )dUi • • • du m = 



I 



r(a + . . . + a m + m) 



til,... ,u n 



Proof of Lemma 1.6.2. Induction on m. 

Note that by this statement we can immediately compute the moments of Fi, 
F 2 , and F 3 . Indeed, the moments of F± is 2<i-dimensional Dirichlet integral. We 
get 

T=i (*)e(a<+a»+i) 
The moments of F 2 are products of 2d Euler beta-integrals, 

<^,n(K:)Wo>-n^ +1M ^ +lw 



Finally, the moments of F 3 are the products of d integrals of monomials (aQ 
over [0,1]. Thus 

(^m^k<) A "))=u A ,, 1 A „ +1 - 

i=l i=l i ^ 
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>\A'+A'' 



The moments of pseudoconvolution are products of moments of factors, see Propo- 
sition 1.4.2. Thus, we obtain (1.18). □ 

Proof of Proposition 1.5.1. We set 

(1.21) I d (a' 1 ,a'{;...;a' d ,a'^)= ^ sgn a ■ J d (a[, a£ (1) ; ...; a' d , a£ (d) ). 

Then (1.18) implies (1.6), and (1.19) implies (1.7). □ 

By the fact that F± and F 2 are invariant under the permutations of {a"}, we can 
write down the following formula for I d , which easily follows from (1.20), (1.21). 

Corollary 1.6.3. For every d = 1, 2, . . . 

(1.22) I d = sgn a ■ F 1 F 2 F 3 CT , 

creS d 

where Fi , F 2 are as above, and 



F^a^a";... ;a' d ,a d ) =5(a[ -a" (1) ,... , a d - a" (d) ) J[ X[o,i] ( a D 

i=l 

7-1 / / // / // \ 



Here 

( , / 1, « e [0, 1] 
X[0 ' 1](w) = \o, «*[0,1] 

is £/ie characteristic function of the segment [0,1]. 
Example 1.6.4. By the formulas above, 

h(a',a") = Ji(a',a") = f 1 Qf 2 Qf 3 , 

where 

/i(a',a") = rWAa'W-Aa") ' ^-2(1 - (a + a")); 
f 2 (a',a") = 4> z -z-i{a')4>-z,z-i{a")\ 
f 3 (a',a") = S(a'-a") X[ o,i ] (a'). 

This statement can be rewritten in the following form. 

h{a', a") = T(t) J ^ z ,{u)^ zl {v) • &_ 2 (1 -u-v) 

u,v 

f ( a' \ , ( a" \ dudvdw 

(1.23) x 



z-z-l y- z z-i ir- 

v uw I V vw I uvw z 



ro£[0,l] 
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Indeed, the value of such distribution on a test function ip(a',a") equals, by defini- 
tion, 



T(t) J <f> z >(u)<l>- Z '(v)-<f>t-2(l-u-v) 

u,v 



X 



H'z,— z- 
m£[0,l] \a',a" 



-i(a')(/)- Z;Z -i (a")i()(a'uw, a"vw)da' da' 



dudvdw 



(1-24) =(/l0/20/3^>, 

see 1.4.2. 

If we want to interpret the integrals in (1.23) over u,v,w as the usual ones and 
not as a formal sign of pairing, we need to impose some conditions; in particular, 
the set {(it, v, w) \ uvw = 0} has to be negligible. From (1.24) we see that if /i is an 
integrable function, then everything works (the expression in brackets is a smooth 
function in u,v,w). Thus, in order to understand the integration in (1.23) as the 
usual one, we may require, for example, — 1 < Rz' < 1, t > 2. We shall discuss 
these problems in more details at the beginning of Chapter 2, where we shall use 
the multidimensional analog of (1.23). 

1.7. Combinatorics of the substantial structures. In this section we give a 
nice combinatorial reformulation of Theorem 1.5.8. 
Let us denote by $ n> d the set of mappings 

¥ >:{l,...,n}-{l,l';...;d,d'} 

subject to the two conditions 

1) <p is injective, i.e. <p(i) ^ <p(j) if i ^ j; 

2) Imp fl {m, m'} ^ for all m = 1, . . . , d. 

It is clear that <& n ,d 7^ if and only if n/2 < d < n. On the set {1, 1'; . . . ; d, d'} 
we have a natural action of the symmetric group Sd' this group permutes the pairs 
(i,i'). The action induces an action of Sd on $ n) d, and the condition 2) implies 
that every orbit of this action consists of exactly \Sd\ = d\ points. 

Proposition 1.7.1. There exists a map 7 n> d of the set & n ,d onto the set of all pairs 
(T, T) where T is a substantial structure with n blocks and d fragments, and T is 
an enveloping structure of T . The inverse image of any pair (T,T) with respect to 
this map consists of exactly one orbit of Sd in $ n ,d- 

Proof. Let us construct 7 n , d. We fix ip G <& n ,d and produce the pair (T, T) = 
7n,d(v ? ) as follows. Let us describe the ith fragment of T and T. 

(1) If (p~ l (i) = then the ith fragment of T has only one (necessarily hook) 
block, and the number of its block equals </? -1 (z'). The enveloping structure T has 
in this fragment one extra linear horizontal block. 

(2) If </?~ 1 (z / ) = then the ith fragment of T also has only one (necessarily hook) 
block, and the number of its block equals </? -1 (i). The enveloping structure T has 
in this fragment one extra linear vertical block. 

(3) If </? -1 (i) > </? -1 (z') then the ith fragment of T has one hook block number 
</? _1 (i / ) and one linear horizontal block number </? -1 (z). The ith fragment of T 
coincides with that of T. 
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(4) If ip~ l (i') > <p~ l {i) then the zth fragment of T has one hook block num- 
ber </? -1 (z) and one linear vertical block number tp~ l (i'). The ith fragment of T 
coincides with that of T. 

Clearly, 7^ is a surjection. Moreover, given a substantial structure T (with 
n blocks and d fragments) together with a numeration of its fragments and its 
enveloping structure T, we can restore <p using the conditions l)-4) uniquely. A 
change of the numeration of d fragments of T exactly corresponds to the action of 
Sd on <& nj d described above. This numeration is not determined by T or T, thus, 
the inverse image of any pair (T, T) is exactly one orbit of Sd- □ 

Let us introduce some notation. Starting from a function (or a distribution) 
F(ri, si; . . . ;rd,Sd) in 2d variables and a map (p G <& n ,d we define the function 
((pF)(xi, . . . ,x n ) in n variables, <p> G $> n ,d, as follows. Let us rename the variable 
ri by Xk if <p(k) = i, sj by xt if <p(k) = j', and let us do this for all ri, Sj such that 
G Innp. Then let us integrate F over all 77, s m such that / ^ Im</?, m' limp. 
The result is denoted by (<pF)(xi, . . . ,x n ). 

For the sake of convenience we also introduce distributions (d = 1, 2, . . . ) 

t d 

(1.25) H d (n, si; . . . ; r d , s d ) = I d (ri, -si; . . . ; r d , -s d ) 

d\ Yl(ri + \si\) 
i=i 

where Id is defined by (1.22), cf (1.13). These distributions are defined in the 

d 

domain where Yi ( r i + \ s i\) 0- Proposition 1.5.1 implies that 

i=l 

d 

(1.26) snppH d c{(r 1 ,s 1 ;...;r d ]Sd)eM. 2d \r l >0, < 0, ^(r, + | Si |) < 1}. 

i=i 

The distributions Id(ri, si; • • • ; r^, s,/) are symmetric under the permutations of 
pairs (r^, (because their moments (1.6) are symmetric). Therefore, by (1.25), 
Hd(ri, si; . . . ; r^, s^) are also symmetric under these permutations. 

Proposition 1.7.2. For any ip G <& n ,d 

(ipH d )(x 1 ,...,x n ) = ^p ( n zz '\T,f) 

where 

(T,f)= ln , d (<p). 

Proof. Follows from the definitions and the fact that Hd(ri, si; . . . ;rd, Sd) is sym- 
metric under the permutations of pairs (ri, Si). 

Finally, combining Theorem 1.5.8 and Propositions 1.7.1, 1.7.2, we obtain the 
main statement of the first chapter. 

Theorem 1.7.3. The correlation functions ph, (xi, . . . , x n ) in the domain where 
all variables are nonzero and pairwise distinct have the form 

n 

Pn Z,) ( x l ■ ■ ■ » x n) = ( l P H d)(xi,... ,X n ). 

d>n/2 n,d 
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2. The Correlation Functions in Positive (Negative) Hyperoctants 



By the words 'positive (negative) hyperoctant' we mean the domain where all 
variables are positive (negative). We have the fundamental relation (see [O, Propo- 
sition 4.6]) 

(2.0) P { : z '\xi, ■ ■ ■ , x n ) = p n - z >- z '\- Xl , . . . , -x n ) 

which implies that we may consider only positive hyperoctants. 

2.1. Why positive hyperoctants are simpler. Remark 1.5.9 shows that the 

support of p n (T) has common points with the positive hyperoctant only if all n 
blocks of the substantial structure T are hook blocks. For each n such structure T is 
unique. Moreover, the only summand of (1.17) which gives a nonzero contribution 
in the positive hyperoctant is such that the enveloping structure T has only extra 
(with respect to T) linear vertical blocks (thus, no extra horizontal blocks). There 
is exactly one such enveloping structure. Therefore, using (1.13), (1.16), and the 
definition of (1.25) we obtain 

Proposition 2.1.1. For x\,... , x n > in the domain where all variables are 
pairwise distinct 

p^ z \xi, . . . ,x n ) = n\ J H n (xi,si;...;x n ,s n )ds 1 ---ds n 

s 1 1 ■ ■ ■ i s n 



(2.1) =t n / ■ rl n (x 1 ,s 1 ;...;x n ,s n )d8 1 ---d8 r 



Si , — fS n 



where the distributions I n were defined in Proposition 1.5.1. 

Thus, we reduced the complicated combinatorial expression for the correlation 
functions obtained before (Theorem 1.7.3) to the integral (2.1), provided that all 
our variables are positive. 

2.2 An integral representation. In this section we shall derive from Proposition 
2.1.1 the following formula. 

Theorem 2.2.1. Let x\, . . . , x n > and Yl7=i x i < 1- Then in the domain where 

UiKjiXi-Xj) ^ 0, 

» n 

p(* x '\x u . . . ,x n ) = t n T(t) J J] </>_ 2 (a,)<Ma, + 1)0-^)0^ + 1) 



CLi,bi 

i=l,... ,n 



i=l 



( 1 \ n 

(2.2) x det (pt-n-i (1 - x ^ + 6 * + !)) II daidhi ' 

\a i + o j + i.; , =l , =l 



Remark 2.2.2. In Section 2.5 we shall prove that the set where n^j^i — x j) = 
is a null set for the correlation functions and that the restriction n^jl^i ~ x j) 
on the domain in Theorem 2.2.1 is irrelevant. 
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Proof. Note, first of all, that the right-hand side of (2.2) is a well-defined generalized 
function in our domain. Indeed, the product 

n n 

Y[(f>-z{a>i)(f>z'{ai + l)<f>-e'(bi)<f> g (bi + 1) • 0t_ n _i(l - ^£;(a; + bi + 1)) 
i=l i=l 

is a well-defined distribution in a^s and b^s, because any point (a, b) belongs to 
at most 2n singular hyperplanes corresponding to the singularities of factors. The 
integral in (2.2) is the value of this distribution on the test function 

det 



ai + bj + 1 



Secondly, the right-hand side of (2.2) depends on z and z' analytically, and so 
does the left-hand side (the correlation function) . That is why in this proof we may 
restrict ourselves to the domain where —1 < Ufc, 3tz' < and t = zz' > n+ 1. Then 
all integrals appearing in this section can be considered as the usual ones, and not 
as formal signs of pairing (cf. the end of 1.6). Moreover, all distributions <p a in 
this proof become just usual integrable functions, and we shall operate with them 
accordingly to this. 

By the explicit formulas (1.20), (1.21) for I d and the definition of the pseudo- 
convolution (section 1.4), we obtain (cf. Example 1.6.4) 

/n n 
Y[ 4>z> (Ui)<f>- Z ' (Vi) ■ (pt-n-l (l - ^2(Ui + Vi)) 
i=l i=l 



Ui ,Vi 



(2.3) x j2 sg^) / n (—) (- ?L -) 

,es n J ti Kuwait)/ 



duidvi 



Wi e[o,i] 

i—l,... ,n 



Let / be a distribution with compact support, and supp / C R+. A function &f 
of one positive variable is called the Stieltjes transform of / (cf. [W]), if 

&f(y) = [ ^-dx. 

J x + y 

X 

Clearly, 

(ef)(ay) = e(f(ay)), a > 0. 
Substituting (2.3) into (2.1) and using this observation, we obtain 

/n n 
11 4> z , ( Ui )<f>- Z , (Vi) ■ (Pt-n-1 (1 - + Vi)) 

i=l i=l 



i 

Ui,Vi 



X 



i=l,... ,n 

(2.4) 

aeS n , 1=1 x \ i ' 

«>ie[o,i] 

i=l,... ,n 

We shall need the following 
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Lemma 2.2.3. 

&<P- z , z - 1 (y) = y- z (l + yy- 1 . 

Proof of Lemma 2.2.3. We shall denote by F(a, b; c; x) the standard Gauss hyper- 
geometric function, 

771/ 7 X r / N f 4>b-l(t)(j)c-b-l(l ~ T) 

F(a,b;c-x) = T(c) j {1 _ XT y dr > 

T 

see [E]. We have 



x + y y 7 1 + x/y 

X X 

= + 1, 1; 1; -1/y) = -(1 + Vy)" 2 " 1 = y~ z (l + y)^ 1 

y y 

where we used well-known following identities 

(2.5) F(a,b;c;x) = F(b, a; c; x); F(a, 6, b; x) = (1 - x)~ a . □ 

By applying the lemma to (2.4) we get 

/n n 
J[(l>z'(Ui)(f>- z >(Vi) ■ 0t- n -l(l " ^(iH + Vi)) 

Ui ,Vi 

i- 

(2.6) 



i=i i=i 



71 \ / \ , s y 1 



0<u>i<l 



duidvidwi 



X 



UiViwf 



We proceed by the following 
Lemma 2.2.4. 

l 



I" l ( ri \( r2 \~ z (l\ r 2 V" 1 dw _ (p- z (u - ri)0 z (r 2 + v) 

J Z ' Z \UwJ \VwJ V VwJ UVW 2 TiV + T2U 







Proof of Lemma 2.2.4- First, let us simplify the expression 

\uw / \vw / V ot/ uvw 
1 / ri y 1 A _ n_\~ 2_1 (^IfLV 2 Ci + IlV" 1 1 



r(z + 1) v-uw/+r(— ^) v uw J + \vwJ v tj^y m>«; 2 

(uw — ri)^ 2-1 (vto + r 2 ) 2_1 

f(=i) r(s + i) 
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Let us make the change of variables w — > y 

uw — T\ (u — ri)y + r\ , u — r\ 1 

y = ; w = ; aw = ay. 

it — r\ u u 

Then we obtain 



/ 



1 

{uw — T\)\ Z (vw + r 2 ) 2_1 . 

aw 



T(-z) T(z + 1) 

o 

1 



/ 



(y(u-r!)) + z 1 ((y(u-r 1 ) + r 1 )v/u + r 2 ) z 1 u - r\ 
T(-z) f(iTl) v~ 



dy 



= «-(« - n)?^±^ / (l + =^Vf' dy 

T{z + 1) J T(-z) \ vr 1 +ur 2 J 
o 

-zt \-z ( vr i + ur 2 ) z ~ 1 1 / u(u-n) 

-u (it — n) + z— rr — —F I 1 — z, —z; 1 — z; 



T(z + 1) r(— z + 1) \ ' iTi + itr 2 

by the definition of the hypergeometric function, see above. Using the identities 
(2.5) we finally get, that our initial integral equals 



_ z (u - n) + z (vn + ur 2 ) z 1 / v(u - ri) V _ (it - n) + z (v + r 2 ) z 1 
U T(-z + 1) I\z + 1) V wi + ur 2 J ~~ T(-z + 1) T(z + 1) vr x + 



+ ur 2 



The proof of Lemma 2.2.4 is complete. 



To complete the proof of Theorem 2.2.1 it only remains to use the result of the 
previous lemma in (2.6) n times taking 

ri = Xi, r 2 = x a -i(i), u = Ui, v = v a -i^), w = wf, i = l,...,n. 

Then we have 

/n n 
II i U i)^-z' (Vi) ■ (pt-n-l (1 - + Vi)) 

i=l i=l 

Ui,Vi 

i=l,... ,n 

1 

x ^2 s S a (^)Y\. < l>-z(ui ~Xi)(f) z (vi +Xi) duidv 



aES n i=l 

We arrive at the formula (2.2) via the change of variables 

di = Ui/Xi - 1 

bi = Vi/xi 



X<j(i)Vi -\- XiUfj^ 



We conclude this section by the following statement which will be used later. 
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Proposition 2.2.5. The right-hand side of the formula (2.2) is an analytic func- 
tion of xi, . . . ,x n in the domain 

n 

{xi, . . . ,x n > 0; < 1}. 

i=i 

Proof. Set 

By the change of variables 



A-i — i-\ x \ a i 



l-i^r* 

we obtain the following formula for the right-hand side of (2.2) 

/i I \\t-\-n — 1 c n n 

rr(t) [ 2 / \{cj>_ z {Ai)<t>-ABi)<t>t-n-i{i - £(A + so) 
[U=i x i j i=1 i=1 

i=l,... ,n 



1 — \x \ I \ 1 — \x 



x det ( -— 3— ^ I TT dAidBi 



z=l 



This formula can be interpreted as the value of the distribution 

n 

0_ 2 (^)0-z'(^)0t-n-l(l-Xl(^ + ^)) 

i=l 

on the test function 

t n r(t) (1 ~ l ^ l)t+ 2 w " 1 ^ ( Ai + T^nl ^ (Bi + ' 



i=1 xf \ l-\x\J V 1-F 

x det 



The test function is an analytic function of the variables Aj, Bi and of the pa- 
rameters xi, . . . ,x n . Clearly, it implies that the value of a distribution that does 
not depend on the parameters on such test function is also an analytic function in 

Xi, . . . , x n . n 

Remark 2.2.6. Another proof of Proposition 2.2.5 can be obtained from Theorem 
2.4.1. 



Proposition 2.2.5 and Theorem 2.2.1 immediately imply 
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Corollary 2.2.7. The correlation functions pn Z \x±, . . . , x n ) are analytic in the 
domain where all variables are of the same sign, pairwise distinct, and the sum of 
their absolute values is less then one. 

In 2.5 we shall see that the restriction n^O^ ~~ x j) can be removed (cf. 
Remark 2.2.2). 

2.3*. Lauricella hypergeometric functions of type B. Let us recall that 
sections marked by * are considered to be optional, see Introduction for details. 

One can notice that the integral representation (2.2) of the correlation functions 
looks like those of multivariate hypergeometric functions. Essentially, the only 
thing that distinguishes the integral (2.2) from a hypergeometric integral (2.8) is 
the determinant 

1 



det 



a% + bj + 1 



Our goal in Sections 2.3 and 2.4 is to get rid of this determinant and to give an ex- 
plicit expression of the correlation functions via so-called Lauricella hypergeometric 
functions of type B. Our reference for multivariate hypergeometric functions is the 
book [AK]. 
Let 

a = (ai, . . . , a m ) G C m , 



m 



&=(&i,... ,& m ) eC 
ceC\{0,-l,...}, 

y = (2/1, ••• ,y m )- 

Then the m-dimensional Lauricella hypergeometric function of type B is defined by 
the series 

(2.7) F^(a,b;c\y) = £ ^hMh \ \ \ M km (brn) km yfl . . . y ,_ 

fci,...,fc ra >0 l C Jfcl + ...+fc m fc l---- fc m. 



rn 



that is absolutely convergent for \yi\ < 1, . . . ,\y m \ < 1, see [AK]. 

This function can be analytically continued into the domain Kyi < 0, . . . , 3?y n < 
by the Euler-Laplace type integral 

(2.8) FH(«,Ms/)=r(c) / n ^ < i - e "<) n ^ 



til,... ,u r 



In the next chapter we shall use another type of integral representations for these 
functions, namely, Mellin-Barnes type integrals. 

Remark 2.3.1. The function for m = 1 coincides with the Gauss hypergeo- 
metric function (cf. [E]), and for m = 2 it coincides with the Appell hypergeometric 
function F 3 , see [AK], [E]. 

Proposition 2.3.2. The function (a, b; c\y) is invariant under the permuta- 
tions of triples (ai, bi, yi) and under the transpositions <-> bi for all % = 1, . . . , n. 

Proof. In the domain \yi\ < 1, . . . , \y m \ < 1 the claim is obvious from the series 
(2.7), and the invariance survives after the analytic continuation. □ 
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Proposition 2.3.3. For all k = 1, . . . , m 



F^\aMc\y)= a -^F^\aMc+l\y) 



dy, 
where 

a = (ai, . . . , a fc _i, a k + 1, a fc+ i, . . . , a m ); 
6 = (6i, . . . , 6fc_i, 6 fc + 1, bk+i, • • • , 6 m ). 

Proof. As in the previous proposition, the statement is obvious in the domain where 
the series (2.7) converges; and the property in question is stable under the analytic 
continuation. □ 

The integral representation (2.2) of the nth correlation function is given by a 2n- 

(zz 1 ) 

dimensional integral. Hence, it is natural to express pn via the 2n-dimensional 
Lauricella function n ' . It turns out to be convenient to introduce certain function 

(zz') 

in 2n variables such that p n will is its restriction to a n-dimensional plane, and 
to express this new function in terms of Fg U \ We shall do this in 2.4, and now we 
introduce another function in 2n variables, which will be extensively used in the 
next section. 

(zz ) 

Consider the function /A (y'i , • • • , y' n ; y'{ , ■ ■ ■ , y'n) in 2n variables defined by the 
following formula 

ft'\y[,...,y' n ,y';,...,y';) 

(2-9) = my }_ yl{) E (- 1 ) Eei n^) 1 " ei ^ / ) eiF l nl ( a£ ' 6£ ; c ^) 



£i=0,l 

i=l,... ,n 



where y[ ^ y'[ for alH = 1, . . . , n; and 

y = (y'i,-- - ,y'n,y",--- ,y'n) 

a e = (1 - ei - z', . . . , 1 - e n - z'\ e\ - z, . . . , e n - z) G C 2n , 
6 £ = (1 - ei - z, . . . , 1 - e n - z; e x - z' , . . . , e n - z') E C 2n , 
c = t — n(z + z' — 1). 

Proposition 2.3.4. The function (j/i, • • • , y' n \ y'{-, • • • , y'n) can be uniquely 

continued to the points where y\ = y" for some i = 1, . . . ,n, so that the result 
will be an analytic function in y[, . . . , y' n ; y'{, ... , y'^. 

(zz') 

The extended function will be also denoted by f n . 

Proof. It suffices to prove that the sum over all in the right-hand side of (2.9) is 
skew-symmetric with respect to the transpositions y[ <-> y'[ for all i = 1, . . . , n. For 
the function F^ n \a e , b e ; c\y), as follows from Proposition 2.3.2, the transposition 
y'i ^ Vi i s equivalent to the transposition of the pairs 

b l) = (l-ei-z',l-6i- z) < — »• (a -z,€i- z) = (a e n+i , b £ n+l ) 

that is equivalent to two transpositions a\ <-> 6|; a n+i <-> 6^ +i (which do not 
change anything, see Proposition 2.3.2) and the change q — > 1 — €j. For the factor 
(yi) 1 ~ ei (y'i') ei the transposition <-> is also equivalent to the change — > 1 — e^. 
But under this change the factor (— 1)^ £i changes its sign. □ 
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Example 2.3.5. Let us consider the case n = 1. Then by applying the L'Hopital 

(zz') 

rule to the definition of /A and using Proposition 2.3.3 we obtain 

f[ zz '\y\y') = F ] §(aMc\y) 

+ y' (if (a', b>- c+l\y)- {1 _ ^ _ ^ («"> &"; c + My)) 

where y = (y', y'), c = (1 — z)(l — z'), and 

a = (1 — z' , —z), b = (1 — z, —z'); 
(2.10) a' = (2-z\ -z), b' = (2-z,-z'); 

a" = (l-z',l-z), b" = (l-z,l -z'). 

2.4*. Expression of the correlation functions via Fg (a, b; c\y). The main 
formulas of this section will play a crucial role in Chapter 4. 

Now we are introducing the function in 2n variables such that the nth corre- 
lation function is its restriction to a n-dimensional plane, see the discussion after 
Proposition 2.3.3. 

Let us define the function n£ z (y[, ... , y' n ; y" , . . . , y") in 2n positive variables 
by the formula 

/n 

CLi,bi 

i=l,... ,n 

(2.11) 

—r — j ) 0t-n-i (l - J] W( a * + X /2) + + 1/2))) J] daidbi. 

J ' i=l i=l 

Comparing (2.11) and (2.2) we see that if n^j^j " x j) 7^ then 

(2.12) R n (xi, . . . , x n ; xi, . . . , x n ) = (xi, . . . , x n ). 
Set c = t — n(z + z' — 1). 

Theorem 2.4.1. For x\, ... , x n > swc/i t/iat ^™ =1 Xj < 1 and n^O^i ~ x j) 7^ 

n 

Pn Z '\ x l, ■■■ , x n) = T(t) Y\4>z-l{Xi)4>z'-l{Xi) ■ c -l(l - |x|) 

i=l 

(2.13) x sgna •/(-')(?/!,... 

7 2/<t(1)? • • • ) Vain))- 

o"6S„ 

where |x| = XlILi anc ^ / or all i = 1, . . . ,n 

1 — lx| 



(2.13') 



Xj 

32 



The function fn ZZ ^ is explicitly expressed via Lauricella hypergeometric functions 
in (2.9), (2.9'). 

About the restriction ri^oO^ ~~ x j) 7^ see Remark 2.2.2. 

Proof. In fact, we shall prove a more general fact. Namely, if Yll^i^'i + x 'l) < 2 
then 

n I 'I -i- I "I 

i#« Vx, • • • y»; •••,<) = r(t) n^-i(^)^'-i«) • ^-i (! - 2 ) 

(2-14) x ^a-fi zz '\y' 1 ,...,y' n ,y^ 1) ,...,y': {n) ) 

oeS n 

where \x'\ = Y^7=i x 'v \ x "\ = Yl7=i x 'li an< ^ for all i = 1, . . . , n 

1 _ |x'| + ]x"| 1 _ |x'| + |x"| 

(2.i4') -h*— , ° 



By (2.12), (2.13) is an obvious corollary of (2.14). 
Observe that (2.11) implies 

F}( ZZ ')( T ' r'-r" 
(2.15) = £ sgna-^K,... ,< (n) ). 

where 

P(^')f X ' rJ.rJI X ») 

/n 
Il0-*(ai)Ma< + W-Abifaibi + 1) 
(2-16) „..,, " 



n n 

x<t>t-n-i (l - £(^(a< + 1/2) + + 1/2))) J] 



dajdbj 



i=i i=i ^ ^ "* 



Oi + &i + 1 



(To obtain this it suffices to develop the determinant in (2.11) and renumber the 
integration variables {b{}.) 

Next we shall get rid of the unpleasant factor Ylii^ + bi + l) -1 (cf. the beginning 
of 2.3). 

Note that (we use the obvious formula <j)' c = <j> c -i) 



n n 1 

k-n-1 (l " £ + V2) + <(^ + 1/2))) = J] 37—^7 

• i t X 4 X 4 

1=1 1 = 1 1 1 

(2-17) x J] - —J (l - J^yMi + V2) + + 1/2))) . 
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Substituting (2.17) into (2.16) and integrating by parts we obtain 

n 1 

R-n* i X 'li ■ ■ ■ > X 'n'i X 'li • •• i X n) = t H T(t) — 

i=l Xi Xi 

rj f d d \ (j)- z {ai)(j) z/ {ai + l)(i)- zl {bi)(j) z {bi + l) 



/I L 
n 
i=i 



da; dbi I di + bi + l 



ai,bi 
i=l,... ,n 



n n 

x (l - ^(xfai + 1/2) + x'Hh + 1/2))) J] datdbi 



i=i i=i 

Simple calculations show that 

d 8\ 0_ z (a;)<Ma; + l)0_^(6i)^(6i + 1) 



dai dbi J a-i + bi + 1 

= <t>- z (a,i)<l> z >(ai + l)(p- z/ (b i )(i) z (b l + 1) 



(ai + 1)6; ai(bi + 1) 



and we arrive at the formula 

n 1 

-^rf Z H X li • • • > X ni X l > • • • ' X n) = t n T(t) / _ „ 

i=l ^ * 

/" n ( z' Z 

X / []0- g (o i )0 g /(Oi+l)0-^(6i)0 g (6i + l)( 7 - — TT - —77 , 



a,i,bi 
i=l,... ,n 



n n 

(2.18) X0 t _! (l - ^(^(a, + 1/2) + x'Ubi + 1/2))) J] da^. 

i=i i=i 

Now we have to remove all the parentheses 



(2.19) 



z' z 



(ai + l)bi diipi + 1) 



in (2.18). We claim that the summand obtained by taking one of two terms in these 
parentheses for alH = 1, . . . , n coincides, up to the factor 

JL z It'I + It"I\ 

r(f) n^-ite^'-iOO • 0c-i (i - 1 'g' ' ), 

i=l 

with the summand of (2.9) where we put e; = if we take the first term in (2.19), 
and 6i = 1 if we take the second term in (2.19). 

From this fact and (2.15) we immediately obtain the main statement (2.14) and, 
consequently, (2.13). 

We shall check our claim only in the case when all e; = 1, i.e., we take only second 
summands in all the parentheses (2.19). For other values of {e;} the considerations 
are quite similar. 
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Thus, we have to deal with the following expression 



(-*) 



IV n IV 

fr wn^ / n^- 2 («o^(^+i)0- 2 '(^)^(^+i)- (] , , 

i=1 X i X i J i=l a i\°i + L ) 

n n 

x t _ x (l - ^(x<( ai + 1/2) + + 1/2))) J] datdbi. 



a,i,bi 
i=l,... ,n 



The change of variables 

«. = (1-^)3 
u * - ^ 2 ; *» 

leads (after some simplifications) to the following expression 



(r(z)r(z')) n f j - xj' ^ c_1 v 



i 1 

n 



\x'\ 


+ 


x" 




2 y 



xT(c) 



/1 t 
H (1 -^) z '(l - UiyH 

i=l 



>/\z-l 



Ui,Vi 

i=l,... ,n 



X t _i (l - y^(ttj + Vi)) Y\ dUidVi 



i=l 



i=l 



where c = t — n(z + z' — 1), and y^, y" are given by (2.14'). 

The last integral is the Euler-Laplace type integral representation (2.8) for 
Fg n \a e , b e ; c|y) where 

y = (y[,-- - ,y' n ;yi,--- ,y'n) 

a e = (-z', . . . , -z'; 1 - z, . . . , 1 - z) e C 2n , 
b e = (-z, . . . , -z; 1 - z', . . . , 1 - z') e C 2n , 



and the sets of parameters coincide with (2.9') for ei = . . . = e n = 1. 
There is only one thing left — to check the following equality 



(T(z)T(z')y 



n 



(- i ) n nAn^i( a; i)^-i«)' 



X; 



Xi 



This can be easily done using (2.14'). □ 

As an application of all this heavy techniques let us consider the case n = 1. 
Then Theorem 2.4.1 and Example 2.3.5 (together with (2.0)) give us the following 
statement. 
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Corollary 2.4.2. Let 

a = (1 — z' , —z), b = (1 — z, —z'); 

a' = (2- z',-z), b' = (2-z, -z); 

a" = (1- z',1- z), b" = (1- z,l- z'); 

c = -t - z - z' + 1 = (1 - z)(l - z'); 
y = (l-l/x,l-l/x). 

For x > we have 

+ y> (Ff (a\ b'- c+l\y)- {1 _ ^ _ ^ (a", b"; c + l\y)) j 

The same formula holds for x < if we replace x by \x\ and (z, z') by (— —z'). 

Comparing this result with Theorem 5.12 from [O] we see that our techniques 
allowed to reduce the order of the multivariate hypergeometric functions involved in 
the expressions for the first correlation function from 3 to 2. Moreover, our formula 
allows to compute the asymptotics of the first correlation function at the origin. 
The result looks as follows. 

Corollary 2.4.3. For x > 

P ( r\*) = ^^+x z - z 'B z Ax) + x z '- z C z Ax), 

where A zz /(x), B zz '(x), C zz i{x) are some functions analytic in a neighbourhood of 
the origin, 

(2.20) A«,(0) = (1^1^11^ > o. 

7rsin7r(2 — z') 

The same formula holds for x < if we replace x by \x\ and (z, z') by (—z, —z'). 

Note that the constant A zz >(0) is invariant with respect to the change of (z, z') 
by (—2, —z'). That means, that the density of the first controlling measure 

o\ \x) = \x\p\ {X) 

is continuous at the origin. 

Sketch of the proof. If x is a small positive number, then y in Corollary 2.4.2 is 

f2l 

large negative. There is a formula for the analytic continuation of F^' = F 3 to 
the domain of large negative values of variables which represents this function as 
a finite sum of expressions of the form 'minus variable raised to some (complex) 
power times a function analytic at the infinity', see [E, 5.11(10)]. Using this formula 
for all three summands in the formula of Corollary 2.4.2, we get the result. □ 

In Chapter 4 we shall consider much more general situation (asymptotics of the 
nth correlation function) , and the n-dimensional variant of the proof above will be 
carried out in all details. 
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2.5. Simplicity of the processes. In this section we shall prove the following 
statement (cf. [O, Proposition 4.2]). 

Theorem 2.5.1. All processes T'zz 1 Q-tc simple. Equivctlently, for any n > 2 the 

set 

{(xi, ... ,x n )e r\ Y[(x l - Xj) = 0} 
is a null set with respect to pn Z . 

Remark 2.5.2. Here is one more reformulation of the theorem (see [O, Proposition 
3.5]). For almost all, with respect to P zz ', points (ce\(3) € O there are no repetitions 
of the type a>i = a i+1 ^ or fa = fa+i ^ 0. 

Remark 2.5.3. Theorem 2.5.1 and Proposition 2.2.5 imply that the restriction 
lli<;j — xj) 7^ in Theorems 2.2.1 and 2.4.1 can be removed (cf. Remark 2.2.2). 

Proof of Theorem 2.5.1. We shall use Proposition 3.5 from [O]. Namely, we shall 
prove that on the diagonal A = {(xi, x 2 ) £ I 2 \ x\ = X2} we have 

(2.22) o- { 2 zz,) \ A = \x 1 \a[ zz '\dx 1 )5(x 2 - xi). 

Proposition 3.5 and Proposition 4.2 of [O] state that this condition is equivalent to 
the statement of the theorem. 

The fundamental relation (2.0) shows that it is sufficient to check (2.21) only on 

A + = {(x!,x 2 ) e I 2 \ xi = x 2 > 0}. 

As was mentioned in Remark 1.2.2 there are three structures with two blocks. 
Namely, we shall denote by Ti the structure with two hook blocks, by T 2 — the 
structure with one hook and one linear horizontal blocks, and by T 3 — the structure 
with one hook and one linear vertical blocks. Then, see Proposition 1.3.2, 

(2.21) a { 2 zz) = a { 2 zz '\T 1 )+ai zz '\T 2 )+a { 2 zz '\Tz). 

We start with T 2 . Let us compute a 2 2 (T 2 ) using Proposition 1.5.2. We shall omit 
the subscript % because our structure contains only one fragment, and % = 1. We 
get 

at'\T 2 )=tC:, a ,,V a a : b [a'I 1 (a\a'')] 

where 

x\ = a, x 2 = b. 
Applying Propositions 1.4.1 and 1.4.3 we obtain 

a 2 zz '\T 2 ) = tC^'a" [a'h(a',a")-6(b-a')} 

= to* [ h (a, y)-^-.6(b-a) + j^rh (b, -a) . 
J a + y \a\+b 

y 

The relation (1.7) implies that the second summand does not contribute to the 
positive quadrant. Moreover, in the first summand we recognize the first controlling 
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measure multiplied by a • 5(b — a) = Xi5(x2 — #1), cf. (2.1). Comparing this to 

( zz '\ (zz'\ 

(2.21) we see, that we have to prove that A + ia a null set for a 2 (Ti)+a 2 (T 3 ). 

( zz') 

Arguing as above, using Proposition 1.5.2, one can show that cr 2 (T 3 ) gives 
no contribution to the positive quadrant. Furthermore, from Proposition 1.5.2 it 

( zz' ) 

is easy to derive that the contribution of a 2 (Ti) to the positive quadrant equals 
(cf- (2.1)) 

2 f h(xi,s 1 ;x 2 ,s 2 ) 

t X X X 2 / 7 ; T7 ; -ds 1 ds 2 . 

J [Xl + s 1 ){x 2 + s 2 ) 

Si ,S2 

In 2.2 (Theorem 2.2.1) we proved that this is the same as the right-hand side of 
(2.2) taken for n = 2, multiplied by xix 2 . Then, by Proposition 2.2.5, it is an 
analytic function in the domain 

{(xi, x 2 ) e I 2 \ Xl, x 2 > 0, Xl + x 2 ^ 1}. 

(Easy to see from (2.1), it is identically zero if xi + x 2 > 1.) That is why the only 

possit 
point 



possibility for a 2 (Ti) to have A + as a non-null set is to have an atom at the 



(x u x 2 ) = (1/2,1/2). 

(zz 1 \ 

If it were the case, then a 2 ' would have an atom at the same point, and 



(zz f ) I (zz') 

a[ (dxi) = / av, (dxi,dx 2 ) 

J X2 



would have an atom at the point xi = 1/2. This conclusion contradicts with 
Corollary 2.2.7. □ 

From Theorems 2.2.1, 2.5.1, and Corollary 2.2.5 we immediately obtain (cf. 
Corollary 2.2.7) 

(zz') 

Corollary 2.5.4. For any n = 1, 2, . . . the correlation functions p n are analytic 
in the domain 

n 

{(xi, . . . ,x n )\ x u . . . ,x n > 0, S ^ j x l < 1}. 

i=l 

Note that all the correlation functions are identically zero as soon as ^ \xi\ > 1. 

Remark 2.5.5. We suppose that by explicit formulas of Chapter 1 a more general 
statement can be proved. Namely, that for all n = 1,2,... the nth correlation 
function is analytic in the domain 



{{x u ... ,x n ) E I n \ ^2\xi\ < 1}. 



i=l 
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3. Lifted Processes 



3.1. Lifting. Lifted correlation functions. Two previous chapters were de- 
voted to the point processes originated from the measures P zz > on the Thoma 
simplex Q as was explained in Introduction (see also [O]). In fact, G. I. Olshanski 
proved, see [O, Theorem 6.1], that all these measures are concentrated on the face 
fi of O. Let us define a covering fi of fio as follows 

Q Q = O x R + 

(we understand this product as the product of Borel spaces). 

There is a natural one-to-one correspondence between fi and the set 

Qf Q = {d x > d 2 > 0, > (3 2 > 0; ^(d, + ^) = s G 

i 

Namely, the map 

f:Q ^%; f:((a\p),s)^(sa\sp) 

is bijective. In the following we identify Qq and Q' Q . 

Starting from any probability measure \i on Qq let us introduce the probability 
measure \x on Qq as follows 

g t-i 

li = [i® Y(t) e ~ Sds ' 

After a while we shall set \i = P zz > , but the general construction is of certain 
independent interest (see, for example, Proposition 3.1.2). 

Following the general scheme described in [O, §4], we associate with /j, a stochastic 
point process on / = [— 1, 1] \ {0}. Let us denote its correlation functions by 

Pn (%\ > • • • ? %n) ■ 

Now we shall construct a new stochastic point process. Consider the set I = 
R \ {0} as the phase space; sets that do not intersect a sufficiently small interval 
(—£, s) as the test sets; finite and countable systems of points in /, such that their 
intersection with any test set is finite, as the configurations. 

To any point (a\j3) £ Qq = O we attach the configuration 

(di,d 2 ,... -/? 2 , •• •)• 

Then our measure p, makes this configuration random. We are interested in the 
correlation functions of the point process P M thus obtained. We shall call the 
lifting of the initial process P M . 

Proposition 3.1.1. All correlation measures p n of 7-*^ take finite values on test 
sets. Furthermore, 

+oo 

^-ye" S p n (£lS"\... ,X n S~ 1 ) — . 
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Under appropriate assumptions about p n , this integral transform can be inverted. 

Proof (due to G. I. Olshanski). Take a test set A of P M and consider a random 
variable Na on Oo defined as follows. The value of Na on (a\f3) e flo is the number 
of ckj's and — fa's which belong to A. By [O, Proposition 4.1] Na is everywhere 
finite, and for 

A = A(e) = [-1, —e) U (e, 1] 

we have 

N A(s) < \ 

For a test A of "P M we define the random variable N A in the same way: it counts 

the number of 6Vs and —/V s in A. 

By definition of correlation measures, the value of the nth correlation measure 
p n on the test set A n is the factorial moment of N A : 

Pn(A n ) = E(N A (N A - 1) ■ • ■ (N A - n + 1)). 

Thus, it suffices to prove, that the usual moments of N A are finite. We may consider 
only test sets of the form 

A = A(e) = (-oo, -e) U (e, +oo). 

For any (a\/3) £ O and s > we have 

N A(e) (sa\s(3) = AT4( es -i)(a|/?) < ^. 

Hence, 

oo 

oo 

<4/£^ e -. ds = % <00 . 

~ e fc 7 T(t) e k 
o 

The formula (3.1) is obvious from the definition of the correlation functions. □ 

In the following sections we shall see that the lifting substantially simplifies the 
formulas for the correlation functions of our processes V zz > . Moreover, it also works 
for the Poisson-Dirichlet processes, see [Ki] for definitions. 

The following claim also follows from section 9.4 in [Ki]. 

Proposition 3.1.2. The lifting of the Poisson-Dirichlet process W(t) is the Pois- 
son process on (0, +oo) with density te~ x fx. 

Proof. Simple calculation. The correlation functions W(t) have the form, see [O, 
Corollary 7.4], 

^_ t n (i-x 1 -...-x n y+ l 

PnK^li ■ ■ ■ i %n) 

•El ' ' ' 
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Then by (3.1) the lifted correlation functions are 



S^ ^ _ s t n (l-^iS ^...-InS" 

p„(xi,... ,x n ) = / ^7— e ds 

T(t) x\---x r 





p g -xi-...-l r 



□ 



X\ ' ' Xj 



~(zz') 

3.2. General structure of p n . From now on we apply the lifting to /j, = P zz '. 
Our main goal in this section is to reformulate Theorem 1.7.3 in terms of lifted 
correlation functions. To do this we need to extend the notion of lifting to the 
distributions with compact supports. In what follows r > is a positive number. 

Proposition 3.2.1. Let . . . , be a distribution in k variables with compact 
support. Then there exists a distribution 

(3.2) 0C T /)(£i,...,&) = Jf(^,...,^?^e-'d8 



not necessarily with a compact support, such that its values on a compactly supported 
test function . . . , £&) is given by the formula 



00 

-T-l 



(3.3) (f^) = J S l^ e - s J /(Ci, ••• ,^6 
G i.- 



<is. 



Remark 3.2.2. Formulas (3.1) and (3.2) represent the same operation if r = t. 

Proof of Proposition 3.2.1. The only thing we need to check is the convergence 
of the integral over s in (3.3). But this becomes obvious if we remember that 
any distribution with compact support is the result of applying some differential 
operator dT 1 • • • d^ k to a continuous function with compact support. □ 

Note that (3.3) is correctly defined if ifj is any bounded test function (not neces- 
sarily with compact support). 

Now we want to extend the notion of pseudoconvolution, see 1.4, to distributions 
which do not necessarily have compact supports. Namely, we shall say that a 
distribution h(£i, . . . , £&) is the pseudoconvolution of distributions . . . , 

and g(£i, ■ ■ ■ , (and write h = f&g) if for any compactly supported test function 
■ ■ ■ , £fc) the following relation holds 



(/i, = J f(£ u . . . , ^ m )g{n U . . . , T]m)lp(tlVl, • • • , £m??m) JJ 



Proposition 3.2.3. For any distributions • • • , £&) one? g(£i, ••-,£&) 

compact supports 

(3.4) C T (f®g) = (£ T f)®g. 

Proof. Obvious. 
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Proposition 3.2.4. The lifted correlation function have the form 

n 

Pn Z '\xi • • • ,X n ) = Yl Yl (<pHd){xi, . . . ,X n ) 
d>n/2 n,a 



where 



H d (n, si; . . . ; r d , s d ) = J — (C t+d I d )(n, -si; ... ;r d , -s d ). 

dlU(Ti + \si\) 



i=i 



Proof. Combining (3.1) and Theorem 1.7.3 we see that it is sufficient to prove the 
following relation 

(3-5) C\<pH d ) = ^^<pH d 

for all <f G $ n ,d- Let us check this for one particular cp, for all others the proof is 
quite similar. Let us take 

such that 



Then, see 1.7, 



(pi%) = i' , i < d 
<p(i) = i — d, i > d. 



t d f 1 

(<pH d ) ( Xl ,...,x n ) = -—j- — — / 

L-_i (Xi+d + \xi\) J 



nr=i (Xi+d + \Xi\) J Y[i= n -d+M + \ X i 

1 n — d + 1 s - - - ■>' d 



d 

x Id(x d +i, -xi, . . . ;x n , -x n - d ;r n - d+ i, -x n _ d+1 ; . . . ;r d , -x d ) |[ dr t . 

i=n—d+l 



Hence, 

/_*— n— 1 



oo 

„i-n-l 





X 



/ - \ 

n™=i (x i+i /s+ \x,/s\) J nsL„-d+i( r > + i a 'i/ s i) 

r n -d+i,... ,rd 

T ( X d+ i -Xi X n -X n - d -X n - d+ i ~ X d\T-T,] J 

x I d , ; . . . ; — , ;r n _d+i, ; ... ;r d , \\dn-ds 

\ s s s s s s J ±J - 



By changing the variables 



Ui = sri, i = n — d + 1, . . . , d, 
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oo 



we arrive at the formula 

{-i) d t d r i 7 s *- d - 1 _ s 

Ui=l(Xi+d + \Xi\) J Ylt=n-d+l( U i + \ X i\) I r ^ 

u n _ d+1 ,... ,u d 

T f Xd+1 -Xl X n -X n - d U n - d+ i -X n - d+1 U d ~X d \ -p-r 

xl d , ; — , ; , ;... ; — , )ds-\\du i 

\ s s s s s s s s J - LJ - 

' i 

that coincides with the right-hand side of (3.5). □ 

It remains to compute (C t+d I d ). 
Lemma 3.2.5. Let a± + . .. + a m = 0. Then 

(*» («,)- ^ (u^ mi (i-f»)= fc - (ui) r(Q ^^^: E " i ° 



Proof. Direct computation. 

oo 

/s a ° v-^ 

4=1 

OO 

o ° i=1 

oo 

= r(a + m + 1) / ^ ao ( s ) eS_Eri1 " i<te ' ^ • ' • ^ m (u m ) 



(mi) • • • 4>a m (u m )e~^^ 



11 i 



□ 



r(ct + m + 1) 

Let us use the explicit formula (1.22) for I d . The previous lemma shows that 

i=l 

Thus, using (3.4), we obtain 

(C t+d I d )(a t 11 a t U---ia , di a'£) = det N{a'^ a'D 

where 
(3.6) 

N(a, b) = (<p z ,(a)<j ) - z ,(b)e- a - b ) (^-^(a)^,^)) (S(a - b) X[ o,i](a)) . 

This gives us, together with Proposition 3.2.4, the main statement of this section. 
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Theorem 3.2.6. 

n 

Pn Z '\xii--- ,X n ) = Yl Yl (<P H d)(xi, . . . ,X n ) 
d>n/2 pe<S> n , d 

where 

~ t d 

H d {r u sf : . . . ;r d ,s d ) = jjj^y. + j^j det N ( r i, Sj); 

i 

and N(a,b) is given by (3.6). 
Remark 3.2.7. Clearly, 

supped C {(n,si;... ;r d ,s d )\ n > 0, s t < 0, % = 1, . . . , d}. 

Remark 3.2.8. Our main achievement is that the distributions I d which are rather 
complicated functions in 2d variables, were 'reduced' to one function N in two 
variables. 

3.3. Lifted correlation functions in positive (negative) hyperoctants. As 

was mentioned before (the beginning of Chapter 2) we may consider only posi- 
tive hyperoctants. By analogy with Proposition 2.1.1 we immediately see, using 
Theorem 3.2.6, that in the positive hyperoctant 

p£ e '\xi, ... ,x n ) = det M(xi,Xj) 



where 



M(.,,)=/-^<fa 

J $ + y 



is the Stieltjes transform of the kernel tN(r, s) with respect to the second argument. 
In fact, we have already computed the kernel M. Namely, we shall derive from 
Theorem 2.2.1 the following statement. 

Theorem 3.3.1. Let xi, . . . , x n > 0. Then 

(3.7) p^ z '\ Xl ,... ,x n ) = det M(xi,Xj) 

where 

rr e -x(t 1 + l/2)-y(t 2 + l/2) 

M(x,y) =t <j>- z {ti)(t>-z'{t2)(t>z'{ti + l)Mt2 + 1) , - , : dhdt 2 . 

J J ti + t 2 + 1 

ti ,ti 

Proof . We shall use the fact, see [KOV], that the set 

n 

{On,... ,x n ) e r\ = 1} 

(zz'\ 

is a null set for the nth correlation function p n . That is why we shall neglect 
this set while making the lifting (3.1). 

Let us apply C f to both sides of (2.2). Using the relation 



!>i_ n _i (l - Y x ^ a i + b i + !)) 



i=i 

we arrive at our assertion. □ 



1 e -E"=i^(«,+bi + l) 



r(f) 
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Remark 3.3.2. An immediate corollary of the last theorem is that the lifted cor- 
relation functions are analytic in the positive hyperoctants. We suppose that using 
Theorem 3.2.6 one can prove that lifted correlation functions are analytic when- 
ever all variables are nonzero (cf Remark 2.5.5). Note that the initial correlation 
functions were not even continuous for some values of z and z' on the borders 
Y^i=i \ x i\ = 1 °f their supports. (For the first correlation function it can be seen 
from the explicit formula in Corollary 2.4.2). 

Remark 3.3.3. Theorem 3.3.1 means that if we restrict ourselves to the behaviour 
of {en} only (i.e. we consider the point process in the restricted phase space M + C I) 
then all correlation functions are given by the determinants of the type (3.7). Such 
processes were considered by several authors, see [Me], [DVJ]. Furthermore, deter- 
minantal formulas for correlation functions appear in some models of mathematical 
physics, see [KBI]. For a more detailed discussion of these processes see [BO]. 

It turns out that the main result of Section 2.4 (Theorem 2.4.1) also has a nice 
reformulation in terms of the lifted processes. It expresses the kernel M(x, y) via 
the Whittaker function W Rifl (x), see [E, chapter 6] for the definition. 

Theorem 3.3.4. Let xi, . . . , x n > 0. Then 

Pn Z '\ x ^ ■ ■ ■ i x n) = det K(Xi,Xj) 

where 

(3.s) 'M- ^^^ i 

ipi(x) = X~^W Z + Z > + 1 , jS_(x), (P2(x) = X~^W Z+Z '-! ; 

2 ' 2 2 ' 2 

and W K ^{x) = is the Whittaker function. 

Clearly, the kernel K(x,y) is real symmetric: 

K(x,y) = K(y,x). 

We call K{x,y) the Whittaker kernel. 
Proof. We shall prove the formula 

(3.9) M(x,y) = {x/y)^K(x,y), 

which by Theorem 3.3.1 proves the assertion. 

The simplest way to prove (3.9) is to apply C t+n to both sides of (2.18). However, 
we prefer to give an independent proof here. 

We start from the integral representation of M{x,y) in Theorem 3.3.1. Note 
that 



(-^ ( " +1/2, -" (,a+1/2, = (^-^) 



e -£c(ti + l/2)-j/(t 2 +l/2)_ 
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Integration by parts gives 



(x - y)M(x, y) = t II e -*(*i+i/2)-v(t a +i/2) 

U(*l)0-*'(*2)^'(*l + 1)^(^2 + 1) 



X 



ii ,*2 

dt[~~dt 2 



h + t 2 + 1 



dtidt-2 



Simple calculation shows that 

9 9\ 0_*(*i)0_*'(* 2 )M*i + !)^(*2 + 1) 



5*1 <%2 J t! + t 2 + l 

= 0- a (*i)0-^(* 2 )^'(*l + 1)^(*2 + 1) 



(*1 + 1)*2 *l(*2 + l) 



Then we get 
(x - y)M(x,y) 

tl t 2 7 



e 2 



r(z)r(z' 



Using the standard integral representation of the Whittaker function (or that of 
the confluent hypergeometric function \l>(a, c; x), see [E, 6.9.4 and 6.5.2]) 

W K ,^(x) = e-*' 2 x» +1 / 2 J ^_ K _ 1/2 {T){l + rY + ^ 2 e-^dT 

T 

we arrive at (3.9). □ 

4. ASYMPTOTICS AT THE ORIGIN 

In this chapter we compute the asymptotics of both lifted and non-lifted corre- 
lation functions when all variables are of the same sign and infinitely small. The 
results can be considered as applications of the explicit formulas obtained in the 
previous chapters. The asymptotics is also of certain independent interest; it gives 
a possibility to consider new (stationary) stochastic processes obtained from the 
points which are 'infinitely close to zero', see [BO] for details. In other words, these 
new processes encode the information about the behaviour of a^, flk as k — > oo. 

The idea of considering the asymptotics of the correlation functions is due to 
G. I. Olshanski. He also computed the asymptotics in the lifted case (Section 4.1), 
and this served as a prompt for our further results in this chapter. 



4.1. Lifted processes. In this section we shall prove the following statement. 
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Theorem 4.1.1. Let . . . , x n > 0. Then 

. , ~(zz')f \ det k(xi/xj) + r(xi, . . . , x n ) 

(4.1) p\ >{x u ... ,x n ) = J - 

where 

(4.2) k(x) = 



sin tcz ■ sin ttz' x 2 — a; 2 



7rsin7r(2; - 2;') - x~^ 

{O (max{xi}) , if z' = z ^ z 

O (max{xi In 2 , i/z' = z£l\Z 

O ^(maxjxi}) 1 ' z 2 ' j , i/m < z, 2/ < m+ 1, m e Z 

as xi, . . . , x n — > 0. 

Remark 4.1.2. The main term of the asymptotics depends only on the ratios 
of the variables. When z 1 = z, the function k(x) is defined as the limit of the 
expression (4.2) as z' — z — > 0. 

Proof. By the determinantal formula for lifted correlation functions proved in The- 
orem 3.3.4, it suffices to compute the asymptotics of the Whittaker kernel K(x,y). 
It turns out that the cases z 7^ z' and z = z' require different treating. The latter 
is the limit case of the former, and the asymptotics for z = z' (Proposition 4.1.4) 
can be formally obtained from that for z 7^ z' (Proposition 4.1.3). 

The claim of our theorem immediately follows from two following statements. 

Proposition 4.1.3. Let x, y > and z 7^ z' . Then 

sin Ttz ■ sin ixz 1 (x/y)~*~ — {x/y) s ~^' r'(x,y) 

K \ x i V) = : 7 7T 1 1=^ 

7rsin7r(z — z') x — y ^Jxy 

1 f sin 7TZ ■ sin tvz' (xly)~^~—{xly)~^~ ., 

— ■ '■ (•'•••'/) 



xy [ 7rsin7r(z - z') (x/y)* - (x/y)~^ 



r'(x,y) 



where 

O (max{x, y}) , if z' = z 7^ z 

O ^(max{i, y}) 1 2 2 ' j , if m < z, z' < m + 1, m e Z 

as x, y — > 

Proposition 4.1.4. Let x, y > and z = z' . Then 

, . sin 2 7T2 lnx — lny r"(x,y) 

K(x, y) = + — -—- 

it 2 x-y y/xy 

sin 2 nz In x — In y 



xy{ tv 2 {x/ y )\ - 



r + r"(x,y)J 



where 



r" 



(x, y) = O (max{i In 2 x, j/ In 2 w}) 
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as x, y — > 

In what follows we shall use the symbol to relate expressions the difference of 
which has the same order as r'(x, y) or r"(x, y) (depending on the case we consider). 

Proof of Proposition 4-1-3. We start with the following formula that follows from 
the definition of the Whittaker function and basic properties of confluent hyperge- 
ometric functions, see [E], 6.5(6) and 6.9(2). For /j, £ |Z 

x- 1 ' V„,„(s) = e-** (jJt^_*(i /2 - K + ft 2f< + 1; x) 



where 

(a) fc 



$(a,c;x) = ^ 



is the confluent hypergeometric function. We need only the first terms of each 
summand of the expression above. We get 

(pi(x) = x^~ Ai{x) +X - 2— Bi(x) 

<p 2 (x)=x s ^A 2 (x)+x^ £ B 2 (x) 

where the functions Ai(x), Bi(x), % = 1, 2, are analytic in a neighbourhood of zero 
and 

M0) = r 4^4, b 2{0) -^-^ 



Then 



r(i-z)' v/ r(i-z') 



r(z)r(z') \^ x- y 

+( } ^ Bi(g)g2(y)-Bi(y)g2(a) 

x-y 

| (a;/y)^-Ai(a;)g2(y) - (x/y)^ A 1 (y)B 2 (x) 

x-y 

(x/y) d ^B 1 (x)A 2 (y) - (x/y) £ ^ L B 1 (y)A 2 (x) 
x-y 



This formula consists of four summands, let us denote them by Si, S 2 , S3, S4 re- 
spectively. 
Note that 

/ \i \ z — z' z — z f 1 _i_ z — z' 1 -j- z — z' 

(xy) 2 x 2 y+ 2 = x 2 2 y 2+ 2 ^ 0. 
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Indeed, if z' = z then 



1 I z — z 1 -j- z — z 

x 2 2 y 2 2 



cc 2 y 2 < max{a;,y} ~ 0, 
and if m < z, z' < m + 1 for a certain m G Z then l 11 ^] < § and 

1 _l_ z — z' 1 -i- z — z' 1 I z — z' 1 I z — z' 1 I*, ~' 

0<x 2 2 y 2 + 2 <x 2 1 2 'y 2 1 2 < (max{x,y}) 1 ^ 0. 
Further, the functions 

A 1 (x)A 2 (y) - A l (y)A 2 (x) B 1 {x)B 2 {y) - B 1 {y)B 2 {x) 

x — y x — y 

are analytic and consequently bounded in a neighbourhood of the origin. Hence, 
(xy) 2 Si ~ and (xy) 2 S 2 ~ 0. 

The expression S3 + S4 can be split into two parts: the first one is obtained by 
replacing in £3 + £4 the functions Ai(-) and Bi(-) by their constant terms Aj(0) 
and -Bi(O), respectively, and the second part is a rest term. 

An easy check shows that the first part is equals to 

(x/y)^ L A 1 (0)B 2 (0) - (x/y)^i4i(0)S 2 (0) 



T(z)T(z')\ x-y 

(x/y)^B 1 (0)A 2 (0) - (x/y)^B 1 (0)A 2 (0) 



+ 



x-y 

sin 7tz ■ sin ixz' (x / y) — (x/y) ~~ 2_ 



7rsin7r(^ — z') x — y 

where we have used the formula 

r(z)r(i - z) = . 

sin 71 z 

As for the rest term, it can be written in the form 

x^y~^a{x, y) — x~' 1 y fJ 'a(y, x) 



R(x,y) = 



x-y 



where \i = s -^- and a(x, y) is an analytic function near (0, 0) such that a(0, 0) = 0; 

its exact form is unessential. We shall prove that (xy)^ R(x,y) « 0. 
Since a(x, y) vanishes at (0, 0) we can write it in the form 

a(x, y) = xa x {x, y) + ya 2 (x, y) 

with certain analytic functions a\{x,y) and a 2 (x,y). With this notation, we have 

S(x,y) 



(xy) 2 R(x,y) 



x-y 
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where 



S(x,y) = x^ +At y^ _At ai(x, y) + x^ +M y^~ M a 2 (x, y) 

-x^ _M y^ _M ai(y, x) - x^ _At y^ +M a 2 (y, x). 

By symmetry, we may assume x > y. Since the function x i— > S(x, y) vanishes 
at x = y, there exists a point £ G [y, x] such that 

S(x,y) S(x,y) - S(y,y) dS(x,y) 



x — y x — y 

The same argument as above shows that 

x 

Note also that 



dx 



Looking at S(x,y) we see that dS(x,y)/dx is the sum of terms each of which is 
equal to x^ ± ^y 2TM multiplied by an analytic function near (0, 0) and possibly also 
multiplied by y/x. Hence, after substitution x = £, each term is ~ 0. □ 

Proof of Proposition 4- 1-4- We shall use the following formula from [E], 6.8(13): 

— x/2 

r 1 / 2 !^ = r (Ml/2 - k, 1; x) In a; 

' r (1/2 -k) V v 7 ' ' y 



x 



+ £ (V2 , ^ [^(V2 - « + r) - 2^(1 + r)]^ 



r=0 



where -i/>(x) = r'(x)/r(x). As in the previous proof, we need only the first terms of 
the summands. We have 

<fii(x) = Ci(x) lnx + -Di(x) 
(p 2 (x) = C 2 (x) lnx + D 2 (x) 

where the functions Ci(x), Di(x), i = 1, 2, are analytic in a neighbourhood of zero 
and 

1 _ -if>(-z) + 2if>(l) 



d(0) 



C 2 (0) 



r(-z) 
-1 



, D 1 (0) = 
, ^2(0) = 



r(i-z) 

Then by definition of the Whittaker kernel 



K{x,y) 



1 



I*(z) 



In x In y 



T(-z) 
-j>{l-z) +2^(1) 

r(i-z) 

x - y 



| D 1 (x)D 2 (y) - Di(y)D 2 (x) 
x-y 

| lnxdQr)!^) -lnyCi(y)£> 2 (a;) 
x-y 

lny£) 1 ( a; )C 2 (y)-ln a; £) 1 (y)C 2 ( a; ) 
x-y 
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Let us denote these four summands by Si, S 2 , £3, S4, and use the same notation as 
in the proof of Proposition 4.1.3. Clearly, 

{xy) 1 l 2 \ hix lny] < max{x In 2 x, y In 2 y} pa 0. 

The functions 

Ci(x)C 2 (y) - Ci(y)C 2 (x) D 1 (x)D 2 (y) - D 1 (y)D 2 (x) 
x — y x — y 

are analytic and, thus, bounded near the origin. Hence, 

(xy)^Si ^ (xy)^S 2 pa 0. 

As in the proof of Proposition 4.1.3, let us split £3 + S4 into two parts. In the first 
one we substitute the functions C(-) and D(-) by their constant terms C(0) and 
-D(O), and the second part is the rest. Using the well-known identity 

ip(w) = ip(l + w) 

w 

one easily checks that 

1 ( lnxd(0)D 2 (0) - lnyCi(0)D 2 (0) ln 2 /D 1 (0)C 2 (0) - lnxD!(0)C7 2 (0) 



F 2 (z) \ x — y x — y 



sin 2 tc z In x — In y 



tv 2 x — y 



Let us denote the rest term by R(x, y). We have to show that 

(xy)^R(x, y) pa 0. 

We can the rest term in the form 

lux ■ b(x, y) — In j/ • b(y, x) 



R(x,y) = 



x-y 



where the function b(x, y) is analytic in a neighbourhood of the origin and 6(0, 0) = 
0. Then we can represent b(x,y) in the form 

b(x, y) = xbi(x, y) + yb 2 (x, y) 

where bi(x,y) and b 2 (x,y) are also analytic near the origin. If we introduce the 
function 

~ 31 13 31 13 

S(x,y) = x 5 j/2 \nxbi(x, y)+x^y^ \nxb 2 (x, y) —x^y? \nyb 2 (y, x) —x^y? \nybi(y, x) 
then we get 

(xy^R(x, y ) = ^l. 

x-y 
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Again, as in the proof of Proposition 4.1.3, we may assume, by symmetry, that 
x > y and apply the mean value theorem. Thus, for some point £ G [y, x] 



(xy)^R(x,y) = ^ y) - 9§ ^ y) 



x — y dx 



Note that y/£ < 1 and 



£ylny^0, | V^lny « 0. 

If we explicitly compute the derivative dS(x,y)/dx at x = £, then we get a sum 
of expressions estimated above multiplied by analytic in a neighbourhood of the 
origin functions. Clearly, this sum 0. □ 

4.2*. Analytic continuation of Fg (a, b; c\y). It turns out that to compute the 
asymptotics of correlation functions in the non-lifted case is much more difficult 
than to do that for the lifted processes. The reason is simple. In the lifted case 
we can use well-known asymptotic formulas for Whittaker functions, while for the 
original processes we have to deal with asymptotics of multidimensional integrals 
- integral representation of Lauricella functions. The asymptotics in both lifted 
and non-lifted cases happens to be the same. One of possible explanations is that 
the asymptotic behaviour really depends only on the ratios of the variables (Re- 
mark 4.1.2), and lifting does not change them. However, the direct proof of this 
coincidence is unknown. 

In this section we present the Mellin-Barnes type integral representation of Lau- 
ricella function of type B and derive certain formulas for analytic continuation of 
F^ l \a, 6; c\y). These expansions will be used for computing the asymptotics. 

The results of this section represent a generalization of the well known formulas 
for the analytic continuation of the Gauss hypergeometric function and Appell 
hypergeo metric function F 3 , see [E], [Mar], [Exl], [Ex2]. 

We start with Mellin-Barnes type integrals. 

Proposition 4.2.1. // a^, hi ^ 0, —1, —2, . . . for all i = 1, . . . , m then 
F% l \a,b;c\y) = ^ r ^ 

n r(a,)r(6 l ) 



+ioo +ioo n r(a, + 8i )T(bi + s i )r(-s i )(-y i ) fli ds 



1 

X 



(2ttz) 



m 

r- I 

m 



- 1 x - ' v 1 ( ' ' s ' 

i=l 



where arg (— yi) < n and the ith path of integration separates the points Sj = 
0, 1, 2, . . . from the points Si = —ai — n, Si = —hi — n (n = 0, 1, . . . ). 

Sketch of the proof. This formula can be found in [Exl, (2.5.6)] and [Ex2, (5.2.3.7)]. 
One can obtain the proof by computing the residues of the integrand at the points 
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Si = 0, 1, 2, . . . for all i = 1, . . . ,m. Then one gets exactly the series (2.7). The 
correctness of this operation can be checked by the general techniques described in 
[Mar]. □ 

( zz') 

Note that if we want to analyze the behaviour of p n at the origin, then, see 
Theorem 2.4.1 and (2.9), we need to know the behaviour of Fg (a, b; c\y) for large 
negative values of variables. The Mellin-Barnes integral representation allows us to 
continue the Lauricella function to this domain. 

Proposition 4.2.2. Let 9fa/i, . . . , 3?y m be negative and sufficiently large; and ai — 
bi $l Z for all i = 1, . . . , m. Then 



F [ s\aMc\y) =^-^ £ II II ^ 

f] F(ai)T(bi) x={i 1 <...<i p }c{i,... ,m} k=1 1=1 

1=1 X={ji<...<j m - P } 



n 



X 



T (K - a i k - ak)T(a ik + a k ) yj T( ajl - b jt - (3i)T(b Jl + ft) 

ai ,a„>0 k=l 



11 ^,Uo,. ^ afe 11 



13!,... ,f3 m - p >0 

1 



(4.1) X 



p m—p 

r(c- EK + a*)- E 

fc=l 2=1 



Remark 4.2.3. This formula represents Fg^(a, 6; c|y) as a finite sum of expres- 
sions of the form 'product of minus variables in some (complex) powers times a 
function analytic at the infinity'. For m = 1 this is the well-known formula for 
analytic continuation of the Gauss hypergeo metric function, see [E], 2.1.4(17): 

F(a,b:c:w) r(b — a) 1 , 
v ' Y 1 = — ^— —F a, 1 - c + a; 1 - 6 + a, w" 1 

r(c) r(6)r(c-a) (-^) a v ; 

+ r f["7 b \,, 1 , b F (b, 1 - c + 6; 1 - a + 6, w;" 1 ). 
r(a)T(c — b) (—w) b 

For m = 2 the formula consists of four summands and is a known example when 
one hypergeometric series (F3 in this case) is continued by another (it will be F2 
here), see [E], 5.11(10). 

Sketch of the proof. The formula can be easily obtained by calculating the residues 
of the integrand from the Mellin-Barnes type integral in Proposition 4.2.1 at the 
points Si = —ai — k, S{ = —bi — k (k = 0, 1, . . . ). The correctness of this operation, 
as in Proposition 4.2.1, is checked by the machinery developed in [Mar]. The set 
X C {1, • • • , m} reflects the fact that for every i e {1, . . . , m} we have two sequences 
of poles, if we use the first sequence (i.e., Sj = — eij — k for some k = 0, 1, . . . ) then 
we put % in X; if we use the second sequence (i.e., Si = —bi — k for some k = 0, 1, . . . ) 
then i el = {1, . . . ,m}\X. □ 

We shall apply the last proposition to the hypergeometric functions used in 
Chapter 2. The restriction — bi Z implies for us that z — z' Z (see (2.9')) 
which is equivalent to z 7^ z' . In case z = z' we have = bi for alH = 1, . . . , m, 
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and all the poles Si = — <Zj — n, (n = 0, 1, . . . ) of the integrand in Proposition 4.2.1 
are of the second order. This 'logarithmic' case can be worked out by the same 
techniques as the ordinary one (i.e., z ^ z'). Below we present the 'logarithmic' 
version of Proposition 4.2.2. Its special case for m = 1 is the known formula [E], 
2.1.4(18). 

Proposition 4.2.4. 

i=i 

1 



x 



m 

V i=i 
where 

( m \ T'(x) 

h ki (a) = 2V>(fc< + 1) - ^(a, + hi) - V> ( c - £(a, + h) J , $(x) = . 

4.3*. Asymptotics of the non-lifted correlation functions. In this section 
we shall prove the following result, cf. Theorem 4.1.1. 

Theorem 4.3.1. Let . . . , x n > 0. TTien 

(z Z "w v det k(xi/xj) +f(xi, . . . ,x n ) 
Pt >{xi,... ,x n ) = J - 

•El ' ' ' %n 

where k(x) is defined in Theorem 4-1-1 and 

{O (max{xi}) , if z' = z ^ z 

O (max{xi In 2 Xi}) , i/z' = z£C\Z 

O ^(maxjxi}) 1 ' z 2 ' j , if m < z, z' < m + 1, m E Z 

as x±, . . . , x n — > 0. 

Remark 4.3.2. The asymptotics coincides with that in the lifted case. 

Proof of Theorem 4-3.1. We shall use the notation from Section 2.4. As we have 
seen (Theorem 2.4.1), 



Pn z '\ x i, •••,<&„) = T(t) Y\4>z-i{xi)4>z>-i{xi) ■ 0c-i(l - \x\) 

i=l 

X ^ Sg n ^-/n 22,) (l/l 5 -- - ,yn,ya(l),--- ,Va(n))- 

where c = t — n(z + z' — 1) and for all i = 1, . . . , n 

1 — \x\ 

Hi = ? \%\ = %1 ~\~ • • • ~\~ %m 

Xi 
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the function / is denned in (2.9). 
Thus, it suffices to show that 

r(f) n<^-i(^'-i«) • 6-1 i 1 - • • • ^ ' • • • ' M 

n 

Y[ k(x' i /x' i ') + f(x' 1 , . . . , x' n ; x'{, . . . , x'n) 

i=l 

(x[ ■ ■ ■ x' n )^ + ± ^ (x'{ ■ ■ ■x'/ l )^ +± ~2 E ~ 



where 



1 _ |s'| + l*"l 1 _ Is'l + lx'' 
y' = -- 2 // = _t 2_ 

t>l „/ ' ill „ll 



X- x\ 



and 



O (max{xj, a^'}) , if z' = z ^ z 

f^J^ x ") = ) O (max{x; In 2 x[, x'( In 2 <}) , if z' = z G C \ Z 

O ^(max{i[, a^i'}) 1 ' 2 * ') ' if m < z,z' < m + 1, meZ 

as all a^, a;" tend to zero. 

We shall give the proof for n = 1, for n > 1 the proof is obtained by applying 
the same arguments coordinate-wise. 

As for the lifted processes, the cases z ^ z' and z = z' are different. First, let 
z ^ z'. Then by (2.9) 

(4.2) f(y',y") = (/4V, 6°; c\y) - y"F§\a\ ft 1 ; c\y)) 

y y 

where c = (1 — z)(l — z'); 

a = (1 -z';-z), a 1 = (-*'; 1 - z); 

b° = (l-z;-z% b 1 = (-z;l - z'); 

y = (y';y")- 

Now we apply Proposition 4.2.2 (or [E], 5.11(10)) to two hypergeometric functions 
in (4.2). We get 

Fi\a", b ";c\y) - ((-v'Y'-H-v")' GlW ' 
H-y'r'-^-y'Y r( r ( 'rjru ~1T G2(1/y '' 

+(-»')'- 1 (-^)' r ^ff_"^ ) G»(l/»'.l/»") 

+(-«y- 1 (-«'0'' r ^"_'ff ( :^' ) G4(i/i/',i/!/-') 
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where G^s are analytic at the origin and Gj(0, 0) = 1 for all % = 1,2, 3, 4. Similarly, 
F l B > , 6 ; c|y) = ^ I (-y ) (-y ) - ^AfTZ^j G x {l/y ,1/y ) 

H-v'Y' i-y'r'- 1 T r^z)r{i % -z) ^ 2(lM 1,v " ] 
H-yr{-yr- x T -^^^ G 3 {i/y\ 1/y") 

where, again, G^s are analytic at the origin and (5^(0,0) = 1 for all i = 1,2,3,4. 
After we plug these expressions into (4.2) and rewrite the result in terms of solely 
x'^s and x"'s, we shall get the expression of the same form (modulo multiplication 
by (x/y)( z ~ z )/ 2 ), as we had for K(x, y) in the proof of Proposition 4.1.3. Then we 
word for word follow this proof, and, thus, prove the assertion for z ^ z' . 

If z = z' then (4.2) still holds, but we need to use different expressions for the 
hypergeometric function, namely, we use Proposition 4.2.4. We get 

Fj?V,&° ;c |y) = r^f-",)^! "]) ( ln ^) H-y")Li(W, W) 

+ln(-y0(2^(l) - 1>{-z) - m)L2(l/y\ 1/y") 
+ ln(-y")(2V>(l) " ^(1 -z)- m)L 3 (l/y\ 1/y") 

+(2^(i) - ip(-z) - mWO-) - ^(i -*)- m)U(i/y\ i/y") 

where Lj's are analytic at the origin, Lj(0, 0) = 1 for all i; and 

FfV^Vly) = Y^r^^ H-y")Lxii/y\ W 

+H-y f )W(l) - 1,(1 -z)- <Mt))L 2 (l/y', 1/y") 
+ ln(-y ,/ )(2^(l) - M-z) - m)Ls(l/y\ 1/y") 

+(2^(i) - - V»(*))(2^(i) - ^(i - *) - m)U(i/y\ i/y"] 

where, as usual L^s are analytic at the origin and their constant terms are equal 
to one. 

Again, after we plug this expressions in (4.2) we shall get an expression which will 
coincide with the expression for K(x, y) obtained in the proof of Proposition 4.1.4. 
We deal with it exactly as we did in that proof, and, finally, prove the theorem. □ 
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